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Sergei Mesenden Schelkunoff 


Sergei Alexander Schelkunoff was born in Russia 
in 1897. In the early years of his primary and sec- 
ondary education he was interested in almost ev- 
erything except mathematics which he would not 
have passed at all if it were not for the kindness of 
his teachers. At the age of 13, however, he came 
across Cauchy’s Algebraic Analysis. In this book, 
Cauchy defined exponential and circular functions 
throughout the entire complex plane as solutions of 
appropriate functional equations, and thus ob- 
tained their essential properties. From then on 
Schelkunoff became an avid reader of the mathe- 
matical writings of the masters. He was profoundly 
impressed by Felix Klein’s Elementary Mathematics 
from the Higher Standpoint. In Weber-Wellstein’s 
mathematical trilogy, he was fascinated with the 
account of the realization of Lobachevskian and 
Riemannian geometries in Euclidean space, since 
one could hardly demonstrate more dramatically 
that these geometries stand or fall together with 
Euclidean geometry. 

Most of his mathematical education was gained 
through such reading, for he had little formal train- 
ing after he left high school. During his early col- 
lege years he spent most of his time in defense 
plants. Then came military training as an officer, 
active military duty, and finally the chaos of revo- 
lution. It was not until late in September, 1921, 
when he landed in Seattle and went to Pullman, 
Wash. that this nomadic period in his life ended. 

By the end of January, 1922, he had learned 
enough English to enter the State College of Wash- 
ington, and in June, 1923, he received the BA and 
MA degrees in mathematics. Then he joined the 


Engineering Department of the Western Electric 
Company and became engaged in experimental 
studies of electromechanical systems. 

In 1926 he was invited to return to the State 
College of Washington. He did and taught mathe- 
matics until June, 1929. In the meantime he re- 
ceived the Ph.D. degree in mathematics from Co- 
lumbia University. 

Dr. Schelkunoff rejoined the Bell Telephone 
Laboratories (an outgrowth of the Engineering 
Department of the Western Electric Company), 
and became a member of the newly organized 
Mathematical Research Department. It was at 
this time that he began his work on electromag- 
netic theory and its applications. He wrote numer- 
ous articles in this field, and is the author of Elec- 
tromagnetic Theory, Applred Mathematics for Sc- 
entists and Engineers, and Advanced Antenna The- 
ory. He is also the co-author, with Harald T. Friis, 
of Antennas—Theory and Practice. 

He is a member of Phi Kappa Phi, the American 
Mathematical Society, and the Mathematical As- 
sociation of America; he is also a Fellow of the 
American Association for the Advancement of Sci- 
ence, the Institute of Radio Engineers, and the 
American Institute of Electrical Engineers. In 
1942 he was awarded the Morris Liebmann Me- 
morial Prize by the IRE for his contributions to 
the theory of radio wave propagation. In 1949 he 
received the Stewart Ballentine Medal for out- 
standing research in communications and recon- 
naissance from the Franklin Institute. At present, 
he is Assistant Director of Mathematical Research 
at the Bell Telephone Laboratories. 
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Microwaves and Mathematics 


Mathematics and mathematicians play an essen- 
tial role in the development of science and engi- 
neering. This is particularly conspicuous in the 
case of electromagnetic theory and microwave en- 
gineering. In the preface to his famous treatise, 
Maxwell refers to Gauss as the man who *... 
brought his powerful intellect to bear on the the- 
ory of magnetism, and on the methods of observing 
it, and [who] not only added greatly to our knowl- 
edge of the theory of attractions, but reconstructed 
the whole of magnetic science as regards the instru- 
ments used, the methods of observation, and the 
calculation of the results... . ” Maxwell also adds: 
“The great success which these eminent men 
[Gauss, Weber, Riemann, J. & C. Neumann, 
Lorentz, etc. | attained in the application of math- 
ematics to electrical phenomena, gives, as is natu- 
ral, additional weight to their theoretical specula- 
tions... . ” But a note of dissatisfaction is sound- 
ed in “There is also a considerable mass of mathe- 
matical memoirs which are of great importance in 
electrical science, but... they are for the most 
part beyond the comprehension of any but pro- 
fessed mathematicians.” 

As a matter of fact, the influence of mathematics 
on the development of science transcends direct 
application of mathematical methods either in es- 
tablishing precise relations between observable 
physical quantities or in solving specific problems. 
Even more important is the role of mathematics in 
the creation of “physical” concepts and suitable 
models to aid our understanding and mastery of 
natural phenomena. First comes a recognition of 
essential similarities in apparently different physi- 
cal phenomena; that is, a discovery of the analogy 
between them. Then follow abstractions and gen- 
eralizations. As mathematical concepts grow and 
age, they gather flesh about them and gradually 
become “physical” concepts. This happened to the 
energy concept and the impedance concept. And 
somewhere creative imagination enters the picture. 

Examples are plentiful. We have Lord Kelvin 
who formulated his theory of electrical discharge of 
a Leyden jar by drawing a parallel between electri- 
cal and mechanical phenomena. We have Faraday 
who created a model of an electromagnetic field 
which is essentially geometric even though he en- 
dowed his “lines of force” and “tubes of force” with 
some mechanical properties. And it takes a mathe- 
matician to make bold use of analogies and pic- 
tures without being ashamed of it. How many 


mathematical equations do we find in Maxwell’s 
Theory and Hertzian Oscillations, a book of over 
one hundred pages by Henri Poincaré? None. 

Microwave theory and engineering are particu- 
larly indebted to mathematics in its various forms: 
to Faraday’s basic mathematical model of an elec- 
tromagnetic field, to Maxwell’s imaginative con- 
cept of displacement current and to his translation 
of Faraday’s geometric language into analytic, 
thus enabling Hertz to calculate the field around a 
dipole and Lord Rayleigh to predict that electro- 
magnetic waves of sufficiently high frequency can 
travel inside cylindrical tubes. While some mathe- 
matical predictions might have been anticipated 
by experimental discoveries, there are instances in 
which this would have been extremely unlikely. To 
establish experimentally that the attenuation of a 
circular electric wave decreases with increasing fre- 
quency would have required a tremendous amount 
of effort and ingenuity and expense to keep the 
wave “clean.” Would anyone have been willing to 
make the investment without some idea of the pos- 
sible outcome? 

Abstract model building in the microwave the- 
ory is not over. About a century ago, Lord Kelvin 
formulated his telegraphist’s equations for cables, 
or electric transmission lines as we now call them. 
Then, in more recent times, the concept of “normal 
transmission mode” was evolved which enabled us 
to think of coupled transmission lines as a certain 
number of uncoupled transmission lines. More re- 
cently still, the concept of waveguide mode was 
further generalized so that we can think of each 
normal mode as the result of superposition of cou- 
pled modes belonging to a rather arbitrarily chosen 
set. In this way it was possible to obtain from Max- 
well’s equations generalized telegraphist’s equa- 
tions for quite general types of waveguides. By this 
generalization, the gap between Kelvin’s theory of 
transmission lines, based on the concept of dis- 
tributed circuit parameters, and Maxwell’s field 
theory was closed. 

The future? Mathematics should continue to ex- 
ert great influence on microwave theory and engi- 
neering. There certainly will be continued problem 
solving. And, who knows, some ingenious mind 
may create a model which will throw a new light on 
diffraction phenomena, or present us with a partic- 
ularly satisfying way of thinking about propaga- 
tion in large waveguides. 

SERGEI A. SCHELKUNOFF 
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Well Done, Ted 


We are all indebted to Ted Saad for his three 
years of tireless work as editor of our magazine. 

It does not take more than a brief glance through 
these pages to see what an exhausting time-con- 
suming job it must be to both edit and coordinate 
the material for such a highly technical and valu- 
able publication. 

Older members will recall that before Ted Saad 
volunteered his services, Al Beck and Bill Mum- 
ford combined their talents to put out previous 
issues. Al served as Chairman of the Paper Pro- 
curement Committee and Bill served as Chairman 
of the Publication Committee. When they were 
forced to resign the positions in 1954 because of the 
pressure of other business, Ted agreed to take over 
both jobs. 

Ted accepted the satisfaction that goes with real 
accomplishment as his sole reward. Our publishing 
budget is presently based on the assumption that 
we will be able to publish this rather expensive 
periodical without paying staff salaries. 

If it were not for members like Ted and Kiyo 
Tomiyasu, who succeeds him with this issue, we 
would indeed be hard pressed. We would be left 
with the unhappy choice of curtailing publication 
or increasing dues to hire a professional staff. 

Obviously, the first alternative is out. To reduce 
publication would be to impede the flow of the life 
blood of our organization. Holding at most one or 
two meetings a year, as we do, we need the 


TRANSACTIONS to keep in step with each other 
during the months in between. 

As for the second alternative, I think we all 
agree that the dues are high enough. 

Just how Ted was able to carry water on both 
shoulders for us for three years and still keep things 
humming at his research center, Sage Laboratories 
in Waltham, Mass., remains a mystery to me. And 
I’m sure Ted’s wife, and their two daughters were 
not exactly crestfallen to learn that Ted had 
handed the reigns to Kiyo and was free to spend 
more time around the Saads’ West Roxbury home. 

An indication of how seriously Ted felt about the 
publishing work is his choice of Kiyo as his succes- 
sor. Ted had decided to step down some time back 
but stayed on until he found a man who shared his 
enthusiasm for the important project. 

So everything has turned out well. Ted, recently 
installed as our Vice- Chairman, has a bit more time 
for outside interests (that is, if he has any besides 
microwave theory and techniques). 

We, in the group, can be thankful that the 
changeover took place so smoothly and that in 
Kiyo Tomiyasu, we have another double-duty man 
who can provide us with the same high quality we 
have grown accustomed to in our publications. 

Ted, from all of us, our lasting thanks and best 
wishes. Kiyo, we only ask that you keep up the 
good work. You certainly have made a fine start. 

Herbert F. Engelmann, Chairman, PGMTT 
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A Modified Equal-Element Band-Pass Filter* 


R. BAWER7 anp G. KEFALAStT 


Summary—A method is presented whereby considerable im- 
provement in the frequency response of a five-stage, equal-element 
waveguide filter can be realized while preserving nearly all the struc- 
tural simplicity of this realization. It is shown that by increasing 
the loaded Q of the center resonant element of a five-stage, equal- 
element filter, the pass band ripple can be appreciably reduced and 
_ the skirt selectivity improved. The modified design also provides a 
simple means of bandwidth adjustment. 


INTRODUCTION 


ABRICATION problems and cost of a waveguide 
band-pass filter are appreciably reduced when the 

- structure is realized as a cascade combination of 1 
identical sections.! Plots of the power transfer function 
show, however, that the pass band insertion loss ripple 
rapidly increases with n. This fact has generally limited 
the use of the equal-element realization to filters with 
less than five sections. 

In the case of a five-stage filter, a considerable im- 
provement can be achieved while still maintaining 
nearly all the desirable features of the equal-element 
realization. It will be shown that by simply increasing 
the loaded Q of the center resonant element, the 2-db 
pass band ripple inherent in the equal-element case can 
be appreciably reduced and the skirt selectivity im- 
proved. This may be seen intuitively by comparing the 
symmetrical low-pass equal-element prototype with that 
of an equal-ripple Tchebycheff response (see Table [). 


DABEE, | 
Response Prototype Ladder Elements 
C le & Lp PAC TCs 
+ db Tchebycheff 1.414 1.318 2.241 1.318 1.414 
Equal-element 1.414 1.414 1.414 1.414 1.414 


It is clear that, except for the value of C3, the five 
identical elements are reasonably close to the corre- 
sponding Tchebycheff values. Thus, one might reason- 
ably expect that upon increasing the center element ca- 
pacity (directly proportional to the loaded Q of the 
band-pass realization), the resulting structure should 
have a reduced pass band ripple and a skirt selectivity 
approaching that of the Tchebycheff unit. 


* Manuscript received by the PGMTT, September 17, 1956. 

+ Emerson Res. Labs., Washington, D. C. Formerly with Melpar, 
Inc., Falls Church, Va. 

t Melpar, Inc., Falls Church, Va. ; ae 

1G. L. Ragan, “Microwave Transmission Circuits,” McGraw- 
Hill Book Co., Inc., New York, N. Y., pp. 681-682; 1948. 


THEORY 


The general expression for the insertion loss of a sym- 
metrical five-section, low-pass filter is given by 


LL. = 14 w?/4[wt(L’?CPC) — w?(2LCy2 + 2LC\C — LC) 
aC eC ar) amet aD 


where 
w=frequency variable in radians per second 
(6 =jb3 Slip 
Cy = Cs, and 
C;=C. 
Since the only effect of multiplying the Z’s and C’s by 
a constant is a change of the frequency scale, choose 
L=(C,=1. Substituting these values into (1), one can 
obtain the frequencies at which the insertion loss ex- 
trema occur: 


w, = 0, 


w= MCE) 2 CT 
SCE) PCa] 


From the above equations it follows that the number 
of insertion loss extrema, hence pass band ripples, will 
depend upon the value of the center element C. 


l| 


Case I. OO) 5 extrema. 
Case II. Ca7), 3 extrema 
(double root at w=1). 
Case III. 2<C<4.07 3 extrema. 
Case IV. C407 2 extrema. 
Case V. C>4.07 1 extremum. 


The theoretical response of a five-section, low-pass 
filter is plotted in Fig. 1 for various values of the center 
element C. To more clearly visualize the effect of this 
parameter, the curves are replotted in Fig. 2. Note the 
rapid reduction in the maximum pass band ripple, the 
broad minimum, and the uniform reduction in band- 
width as the capacity of the center element is increased. 

The skirt selectivity of a modified equal-element filter 
has been compared to that of a Tchebycheff filter having 
a comparable pass band ripple. Apart from a slightly 
lower insertion loss for w>>1, the general behavior of the 
modified filter, for losses less than 20 db, follows the 
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Insertion Loss in db 


nits 


Frequency in Radians/Sec. 


Fig. 1—Theoretical response of a five-section, low-pass filter for 
various values of the center element. 


Max. Passband Insertion Loss Ripple In db 


3-db Bandwidth In MC/ Sec. 


Center Element Capacity In Farads 


Fig. 2—Maximum pass band insertion loss ripple and a 3-db band- 
width of five-section, low-pass filter for various values of the 
center element. 
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Relative Insertion Loss in db 


5.1 52 5:3 5.4 5.5 5.6 5.7 5.8 59 6.0 
Frequency in kmc/sec 


Fig. 3—Experimental response curves of a five-section, equal-ele- 
ment filter illustrating the effect of increasing the center resonant 
loaded Q. 


Tchebycheff response so closely that no useful purpose 
is served by reproducing the curve. 

A typical response of a modified five-section, equal- 
element band-pass filter is shown in Fig. 3. The im- 
proved response was obtained by simply replacing the 
inductive posts forming the center resonant element 
with posts of larger diameter. 


CONCLUSION 


In addition to retaining nearly all the simplicity of 
the equal-element realization, the modified design af- 
fords a convenient means for simply, accurately, and 
systematically varying the bandwidth with negligible 
pass band distortion. As seen in Fig. 2, the bandwidth 
may be varied as much as +10 per cent while still main- 
taining a ripple of less than 4 db. 

Although this technique has not been applied to 
higher order equal-element filters, it is felt that the gen- 
eral approach is applicable. 


COS Swe 
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Application of Rayleigh-Ritz Method to Dielectric 
Steps in Waveguides* 


R. E. COLLINY anp R. M. VAILLANCOURTT 


Summary—The Rayleigh-Ritz method is applied to obtain ap- 
proximations to the first NV eigenfunctions and corresponding eigen- 
_ values in an inhomogeneously filled rectangular waveguide. These 
approximate eigenfunctions are then used to obtain a solution for the 
reflection and transmission coefficients at the junction of an empty 
_and partially filled waveguide. Theoretical and experimental results 
are given for a dielectric slab which extends completely across the 
broad dimension of the guide, but only partially across the narrow di- 
mension. The experimental values are within the experimental error 
of the computed values obtained by considering the dominant mode 
and only two evanescent modes. 


INTRODUCTION 


AVE propagation in a waveguide inhomogene- 

\ \ ously filled with a dielectric has been studied by 
many authors.! As a general rule, the modes are 

more complex and transcendental equations have to be 
solved, in order to find the propagation constants of the 
various modes. This has led several authors to consider 
the application of variational methods for obtaining ap- 
proximations to the eigenvalues.?-4 By means of the 
Rayleigh-Ritz method (hereafter called the R-R meth- 
od), one may obtain approximations for the frst NV ei- 
genvalues, and also for the first V eigenfunctions.®® In 
this paper, the discontinuity between an empty and an 
inhomogeneously filled rectangular guide will be studied 
using the R-R method. The following procedure is used, 


1) The first V eigenfunctions in the inhomogeneously 
filled guide are approximated by the Rayleigh- 
Ritz method. 

2) Equations expressing the continuity of the tan- 
gential field components at the junction are then 
easily written down and solved for the reflection 
and transmission coefficients. 


* Manuscript received by the PGMTT, October 16, 1956. 
+ Canadian Armament Res. and Dev. Establishment, Valcartier, 


1L. G. Chambers, “Propagation in waveguides filled longitudi- 
nally with two or more dielectrics,” Brit. J. Appl. Phys., vol. 4, pp. 
39-45; February, 1953. (This is a review article containing sixteen 
references. ) 

2 L.G. Chambers, “Compilation of the propagation constants of an 
inhomogeneously filled waveguide,” Brit. J. Appl. Phys., vol. 3, pp. 
19-21; January, 1952. 

3 L. G. Chambers, “An approximate method for the calculation of 
propagation constants for inhomogeneously filled waveguides,” 
Quart. J. Mech. and Appl. Math., vol. 7, pt. 3, pp. 299-316; Septem- 
ber, 1954. 

4 A.D. Berk, “Variational principles for electromagnetic resona- 
tors and waveguides,” IRE Trans., vol. AP-4, pp. 104-111; April, 
1956. 

5 R. Courant and D. Hilbert, “Methods of Mathematical Phys- 
ics,” Interscience Publishing Co., New York, N. Y., 1st English ed., 
ip: 175% 1953. 

6 R. Weinstock, “Calculus of Variations,” McGraw-Hill Book 
Co., Inc., New York, N. Y., 1st ed., chs. 7-9; 1952. 


In a rectangular guide inhomogeneously filled with a 
dielectric slab, as in Fig. 1, the two sets of fundamental 
modes are the longitudinal section electric and magnetic 
modes (LSE and LSM modes), having the electric and 
magnetic vector, respectively, contained entirely within 
a longitudinal section. In an empty guide, the TE and 
TM modes may be derived from a magnetic and an 
electric Hertzian potential having only a longitudinal 
component respectively.? By analogy with this problem, 
it is readily seen that the LSE and LSM modes may be 
derived from a magnetic and an electric Hertzian po- 
tential, respectively and having a single component di- 
rected normal to the dielectric-empty guide interface. 
These two sets of modes form a complete set in which 
any arbitrary field distribution may be expanded.® 


7 ¢ 


zx x 


(a) (b) 


Fig. 1—Inhomogeneously filled rectangular waveguide. 


For a general inhomogeneously filled cylindrical 
guide, the division into LSE and LSM modes is not 
possible. 

When an Ao mode is incident at the junction of an 
empty rectangular guide and a guide partially filled, as 
in Fig. 1(a), both LSE and LSM modes are excited. In 
the case of Fig. 1(b) with an Ayo mode incident, only H,, 
modes are excited. The case of Fig. 1(a) is considerably 
more complex and, therefore, was chosen as a good ex- 
ample to which the R-R method could be applied. 

This same type of dielectric step discontiuity has 
been treated in a recent paper by Angulo, using the cor- 
rect expressions for the eigenfunctions and a variational 
method for evaluation of the equivalent circuit param- 
eters.? However, the coupling of the LSE modes by the 


7J. A. Stratton, “Electromagnetic Theory,” McGraw-Hill Book 
Co., Inc., New York, N. Y., sec. 6.1; 1941. 

8 J. Van Bladel, “Field expandibility in normal modes for a multi- 
layered rectangular or circular waveguide,” J. Franklin Inst., vol. 
253, pp. 313-321; April, 1952. 

9C. M. Angulo, “Discontinuities in a rectangular waveguide par- 
tially filled with dielectric,” IRE Trans., vol. MTT-5, pp. 68-74; 
January, 1957. 
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step when a LSM mode is incident is neglected. There- 
are only approximate and 
e. (tb or t=), where f 
is the slab thickness and 6 the waveguide height) when 
the amplitudes of the coupled LSE modes are small. 

Collin 
f the 


fore, the results obtained 


valid for small discontinuities; 7 


This same step has been analyzed by R. E. 
(unpublished work) as well as the simpler case of 
H-plane dielectric step.!° 

Use of the R-R method has the great advantage of 
avoiding evaluation of many complicated expressions, 
especially for cases where a waveguide cross section is 
divided into more than two regions by dielectric media 
of different dielectric constant. 


WAVE EQUATION FOR THE MODES 


In order to bring the problem being studied into the 
class that can be handled by the normal Sturm-Liouville 
theory, the rectangular guide will be considered as filled 
with a lossless dielectric material with a dielectric con- 
stant which is a continuous function of the coordinate y, 
but independent of x and z. From the theory for the 
continuous case, one may readily pass to the case where 
the dielectric constant is a discontinuous function of y; 
e.g., a slab filling. With a varying dielectric constant, the 
wave equation satisfied by the Hertzian potentials are 
modified and therefore will be derived here. 


LSM Modes 


Maxwell’s equations in a source free medium are 


— = 
VX A = jwenk, (1a) 
=> — 
VX E= — jouH, (1b) 
— 
V7B2— 0; (1c) 
Vien: (1d) 


where x is the relative dielectric constant and here is 
considered as a function of y. 
By virtue of (1c), one may take 


— 


—> 
H = jweV X Us (2) 
= . . 
where IIz is an electric Hertzian potential with a y com- 
ponent only. From (1b) 
= — 
Wipe — eV 


which integrates to 


=> _s 
E = ko’llg + Vo. (3) 
From (1a) 


> — = => 
VX V X Ile = VW Ie — Vile = kez + «Vo 
—_ 
= kko"z + Vib — OVk. (4) 
10 R, E. Collin and J. Brown, “The calculation of the equivalent 


circuit of an axially unsymmetrical waveguide junction,” Proc. IEE, 
vol. 103, pt. C, pp. 121-128; March, 1956. 
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> 
As yet Vo and V-IIg are unspecified and one may, 


therefore, put 


— 


VV-Ilz = Vk 
which apart from an irrelevant constant integrates to 
—- 
V- Ie = xd. (5) 
ae . 
From (5), ¢@=x~!V-IIg and hence the wave equation for 


= . 
IIx, t.e., (4) becomes 


— — = 
V’Ilz + KRo*l pg = K 'VkV - Ip =a) (6) 
while (3) for the electric field becomes 


— 


— — 

E fo Ue + V(k-1V- Tx) 

Vv: TeV! + eNO sit — eee Ty 
KRo” 2p on wav: Te — KV: Tnx 


=x 'VXVX ae (7) 


this latter result following from the wave equation. 
From (7), one sees at once that (1d) is satisfied, since 


4 —> 
V«eE=V‘(VXV Xz) =0. 


For a variation with x according to sin mx/a and ex- 
ponential z dependence, the solutions to (6) are of the 
form 

—* ea TX 


Ig = ty sin — Pp(y)e** (8) 
a 


where Wz is a solution of the Sturm-Liouville equation 


ar dk dwn 1 
a kg a 
a 2 at (et +7) ve (9) 


or equivalently 
d 1 dbz 
ae 
dy « dy 


1 es 
(«te Spee st hi 7 \e =) 
K a” 


Eq. (9) has an infinite number of solutions Wz, with cor- 
responding eigenvalues y,?. These solutions form an 
orthogonal set with respect to the weighting function 
x! and may be normalized so that 


b 
| VanWask ‘dy = Ons (10) 
0 


where 6,, is the Kronecher delta and is equal to unity, 
if m=s and zero otherwise. When xk is a continuous func- 
tion, both Wz, and dWz,/dy are continuous. 

When «x is the discontinuous function, 


KY) =o = Ko Uy — 2) 
where U(y—14) is the step function 


(11) 


OPA, 


12 
RAs ras (12) 


uy = 4 
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Wr, and «-!(dWr,/dy) are continuous in order that the 
tangential field components should be continuous at the 
interface. At y=o, b, dWz,/dy vanishes. The term 


in the differential equation (9) becomes 


dbp 
K (1 — ko _ 
( iy t) 


) 


where 6(y—¢) is the Dirac impulse function. The term 
dz/dy is discontinuous at y=t and hence, the second 
derivative also has an impulse discontinuity at y=. The 
physical reason for these discontinuities is the polariza- 
tion charge in the dielectric. 


LSE Modes 


The longitudinal section electric modes do not have a 
component of electric field parallel to V« and, therefore, 


= aay es 
V-«E = x«kV-E4+ E-V« = 0 


>! 
gives V:-E=0. For these modes, one may take 


> 


: a 
E = — jouwV X Ila (13) 


—S 
where I] is a magnetic Hertzian potential with only a 
y component. From (1a), 


= as => 
V X A = kk X lly = ko2V X «ily, 
since 
= = =53 = 
W S< Pl lye = PAY OX AG SOA AY DI iia, 


This equation integrates to 
= + 
ia KRo7I ag + Vv. 
From (1b), 
> — — 
VV -Ilu — VIL ar = KRo?l1 az + Vv. 


— => 
Let VV - Il. = Vv and the wave equation for II becomes 


— = 
Villa + xko*llu = 0, (14) 
= ty 
while the equation for H becomes 
— as = => , 
(ef = KRo*IL ar 4. VV - Hn = WSCA DK Il y. (15) 


With an x variation according to cos tx/a and exponen- 
tial z variation, the solution to (14) is of the form 


——» ——> eae 
Ily = i, cos — Wu(y)etF, (16) 
a 
where Wy is a solution of 
du 2 ‘ 
poset (tot =" + 8°) ba = 0 (17) 
dy? a 
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Eq. (17) has an infinite number of solutions, Warn, with 
corresponding eigenvalues, 6,?, and these solutions can 
be chosen to form an orthonormal set such that 
b 
if VW usdy = Ons. (18) 
0 


Both Wa, and dWy,,/dy are continuous, irrespective of 
whether x is continuous or not. At y=0, 6, Warn vanishes. 


MINIMUM CHARACTERIZATION OF THE EIGENVALUES 


LSM Modes 


If (9) is multiplied by x Wz and the term involving 
the second derivative integrated by parts once (this 
term vanishes, since zg and x—!(dWz/dy) are both con- 
tinuous and dWz/dy vanishes at y=0, 0), one gets 


b b d 2 
rf kK Wn2dy Sa f (=) 
0 0 dy 


Eq. (19) is a variational expression for the propagation 
constant y”. An extremisation of this equation by that 
class of functions ¢(y), which are continuous with at 
least a piecewise continuous derivative and orthogonal 
to the first K—1 correct eigenfunctions Wx;, with re- 
spect to the weight factor x~!, yields an upper bound on 
the true eigenvalue yx”. By means of the expansion 
theorem for a complete set of functions, one may write 


oO 


Qo = »y, OnWEn- (20) 

0 
The orthogonalization conditions give 
an = 0; n=0,1,°-: K —1. (21) 


Substituting into (19), gives 


(2) b ee) a) ad i OG i 
Poets for yi 
0 


K s=K n=K dy dy 
T 
a KRo? _—. :) WasWin a,a,dy 
a 
00 20 d a b 
= yy > asAn {(-¥e. Me ) | 
s=K n=K SPL AD 


2 awa 
— KO Wis ——— 
if | # dy 
TT 
=P («i re “) Yb | ay : 
a 


The integrated term vanishes and, using (9) and (10), 
the result is 


(22) 


n=K 


ie) ie} 
of Deda? = Dea yee. 
n=K 


180 


Assuming that the eigenfunctions Wz, have been ordered 
so that yo?<yi?<y22 + + + <yx?, the result (22) may be 
written as 


> an” An” : 


a) 


ae Gn *¥n* os YK’) 


* n=K - n=K+1 uy 
sf yg et (29) 
> Gas poe Tee 
n=K n=K 
since ¥,2>Y¥x? for n>K. Only when ¢=Werx, will 


yy? =x’. In general, the approximate eigenvalue is too 
large. A suitable series of functions to use for this ex- 
tremisation are the corresponding eigenfunctions for the 
empty guide. For the LSM modes these are 


orn = ee GOS e— aie n= (); ib SOR 24. 
r ie : 4) 


where €5, is the Neumann factor, 


‘. n = 0 
€on = 
Us ph SY, 
LSE Modes 
For the LSE modes, the variational expression cor- 
responding to (19) is 


»f hae i fe oad) 
0 0 dy 
Pe («te 2 “) vac} dy = 0. 
a? 


As for the previous case, the approximation to the Kth 
eigenvalue by that class of functions which are continu- 
ous and vanish at y=0, 6, and are orthogonal to the 
first K —1 eigenfunctions Wy; is from above. A suitable 
set of functions for this extremisation are again the cor- 
responding functions for the empty guide, 2.e., 


os tes 
oun = eo: iil tli Meta 


The above variational expressions are also valid when 
k is a discontinuous function of y. 


(25) 


(26) 


THE APPROXIMATE EIGENFUNCTIONS 


This section will consider the solution for the first 
N-+1 approximate eigenvalues and corresponding ap- 
proximate eigenfunctions for the case of the LSM modes. 
In the previous section, it was shown that the extremisa- 
tion of (19), with respect to functions which were 
orthogonal to the first A—1 true eigenfunctions, gave 
an upper bound on the Ath eigenvalue. Since one does 
not know the true eigenfunctions, the class of functions 
to be used for the Ath extremisation will be made 
orthogonal to the first K—1 approximate eigenfunc- 
tions. It may be shown that this procedure also yields 
an upper bound on the Ath eigenvalue.!! The proof is 


11 Courant and Hilbert, op. cit., ch. 6. 
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based essentially on the principle that the class of func- 
tions which are used for the first NV+1 approximate 
eigenfunctions is a narrower class of functions than the 
complete set of true eigenfunctions. 

The substitution of a series of the functions dz, into 
(19) and subsequent extremisation leads to a matrix 
eigenvalue problem. The resultant matrix is a symmetri- 
cal real matrix whose eigenvalues are approximations 
from above to the first N+1 eigenvalues. For each 
eigenvalue, a solution for an eigenvector exists and the. 
totality of eigenvectors obtained form an orthogonal set 
with respect to suitable weighting factors. This latter 
result follows from the well-known theory of real sym- 
metrical matrices.!? For this reason, the orthogonaliza- 
tion conditions, which were originally imposed upon the 
functions ¢z,, may be dispensed with. 

From this point on, Wz, andy,” will be used to denote 
the nth approximate eigenfunction and eigenvalue, re- 
spectively. For the Kth approximate eigenfunction take 


N 


Wax = >» OnKPEn 


n=0 


(27) 


where a@,, are unknown coefficients to be determined 
subject to the normalization condition 


b N N 
it x 'Wax*dy = MS, ~ GnkasKP sn = 1 (28) 
0 s=0 n=0 
where 
b 
Pos = Pa =? it Kk benobrs@dy. 
0 
Substituting into (19) gives 
Ag : dors dorn 
bel 
n=0 ead 0 dy dy 
PS: 
— ( xko? — a “ae vs) pushEn¢ a 
a 
= stationary quantity. (29) 
Let 
b 
i oe fades dgzn 
0 dy dy 
2 
= («te = )onebent dy a ee = Lae (30) 
a 
Thus, 


DD asxanx(Tsn — YK?Pen) = stationary quantity. (31) 


s=0 n=0 


Equating all 0/da,, equal to zero for n=0, 1, - IVa 
yields the following set of homogeneous equations 


N 
De Gal ep oe VR Lin) = 0; 


n=0 


s='0,1,---¥. Gy 


2 C. G. Montgomery, R. H. Dicke, and E. M. Purcell, “Princi- 
ples of Microwave Circuits,” M.I.T. Rad. Lab. Ser., McGraw-Hill 
Book Co., Inc., New York, N. Y., vol. 8, pp. 405-409: 1948. 
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For a solution, the determinant must vanish and this 
results in N+1 roots for yx? which are the N+1 ap- 
proximate eigenvalues. For each root, say Yx?, a solu- 
tion for a@,« can be obtained, this solution is unique 
when subjected to the normalization conditions (28). 
The set of coefficients a,« are orthogonal to the set arr 


with respect to the weighting factors P,, for KAR, 1.e., 


N WN 
Ss >) OnRAsKP sn = ORK. 


n=0 s=0 


(33) 


The proof is given in the Appendix. 
For the LSE modes, the set of homogeneous equations 

obtained are 

ib. he dace Lil 


SZ, bnk(Qsn aa BK76sn) = 0, (34) 


n—1 


where 


Orn a ONS ag f {ae ate 
o (dy dy 


a) 
= («te aa :) baba} dy. (35) 
a? 


The vanishing of the determinant yields the first NV ap- 
proximate eigenvalues. The corresponding eigenvectors 
define the corresponding approximate eigenfunctions. 
The coefficients 6,, are subjected to the normalization 
condition 


N 
DY bnx?= A, 1 anna PAS En I 
n=1 

and satisfy the orthogonality conditions 


N 
Sy brKonr aa 0, 


n=1 


z(t" - “Jat x. coe 4, (te! -*) 
a~ 


n=1 


RAK. 


Jweol'o(1 —= R) ie Ay a JME ys Ary 4/— cos” COS geen ay = 


‘At 2 nm N ni ie nw 
cay) GE, (/~ Sint ve —— A, > Silvey) 
a b b x ab b b 
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MATCHING OF FIELDS AT JUNCTION 


Consider the junction of an empty and inhomogene- 
ously filled rectangular guide as illustrated in Fig. 2. 


t 
A 
ta 
a ~~ A, TINO 
b t RA, «sw 
eta 25 s 
foe | poste) 


Fig. 2—Junction of an empty and partially filled rectangular guide. 


Let an Hy mode be incident from the empty guide. The 
higher order modes excited will consist of an infinite 
number of the LSM and LSE modes. It will be assumed 
that only the dominant mode, 7.e., the first LSM mode, 
propagates in either the empty or partially filled guide. 
For an approximate solution, only a finite number of 
modes are considered and in the partially filled guide 
these will be taken as the approximate eigenfunctions 
as obtained by the R-R method. At the junction z=0, 
the tangential field components are made continuous 
and this leads to four simultaneous equations which 
must be solved for the reflection and transmission co- 
efficients. The fields are obtained from the two Hertzian 


— — 
potentials Ilz and Ily by means of (2), (7), (13), and 
(15). There are two different expressions for /, and H, 
which, however, yield the same results. In writing the 
continuity equations for the transverse field compo- 
nents, the functions of x and other common factors will 
be deleted to save space. The following equations, ex- 
pressing continuity of #,, H,, E,, and H,, respectively, 
are obtained 


€0n nr 


Ne 
= Br (vx? = ; Se Onkk * aC OSs rye (36) 
K=0 a7 n=0 b b 
No 672 Oy $1 
~~ —G COS ny 
st ab b b 
NW’ Dr ST 
= Jwe€ > Bryx« = OnK == COS oy = Ss ane y= bsK COS yaa): (37) 
n=0 Ken, UE s=1 b b 
nT De nt 
= Jo » DrBr s bak — sin oo = 2 Br = = QnKK — sin — Vv, (38) 
(Je n=1 ab b b ‘ 
N By ni 
eas Jae ae = “5 2 bnk \/- Ah ees (39) 
iil) n=1 b b 
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where R is the complex reflection coefficient, An, Bn, Cn, 
and D, are unknown amplitude coefficients, and 


ie . 
zi =e 5 = ko’. 
a b? 
The amplitude of the dominant mode on the output side 
of the junction, 7.e., for z>0, is Bo. 
Multiplying (36) and (37) by 
/eon nr 
—= cos == 
b 


b 


and (38) and (39) by 


Ve ae eae ie: 


for. 1=0; 1, , N, in turn, and integrating from y=0 
to y=b, converts these equations to the following 
algebraic equations 


N 
—koX(1 + R)Ao = >) Bx 


K=0 
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AN EXAMPLE 


The reflection and transmission coefficients were cal- 
culated for an Hy) mode incident on a partially filled 
guide, as illustrated in Fig. 2, for values of t/b ranging 
from 0 to 1. The free space wavelength was 3.14 cm, the 
dielectric constant was 2.52, and the internal waveguide 
dimensions were 0.9 X0.4 inch. The dominant mode and 
two evanescent modes (one LSE and one LSM mode) 
were taken into account in both the empty and par- 
tially filled guide. This led to a sixth-order determinant 
which, however, had fifteen of its elements equal to 
zero, and was readily reduced to a third-order determi- 
nant. This latter determinant gave (1—R)/(1+R8) as 
the ratio of two second-order determinants. In Fig. 3 
opposite, the modulus and phase angle of the reflection 
coefficient are plotted, while the phase angle of the 
transmitted wave is plotted in Fig. 4. The computed 
values of the transmission coefficient phase angle are not 
very accurate, because of their small absolute value. The 


DS AnkP no, 


n=(0 


. (40) 
Tv 
(De Ses =)4, as; Br (vx: Fe “\> > askP. sny part ) a . N, 
a K=0 a*/ 520 
N 
Tod — R)Ay = >, Brvedor, 
K=0 
(41) 
i nT , aM nr XN 
JHE. nA n sales ea G, = — JED Sy: BRYKGnk ar by ea Drbnk, Ny Ls 2 ‘ N, 
ab K=0 a0 K=3 } 
nn” site ie & N 
PONG nly San Gees etal JOU sy DrBrbnx — aes >t Br DS: Om d on 
ab K=1 a0 K=0 i 
where (42) 
2 Ae 
ho. fhe eS ko6 8n: 73 ral ) alee ip hy : N, | 
n* ar? a~ J 
2 N 1 
Ga ite =) oo yD (6x? FR =) Onk, I= i; 2 “Vr wn : (43) 
a K=1 a? 


This latter system of equations may be written in 
matrix form as a set of four homogeneous matrix equa- 
tions. For a solution for the amplitude coefficients, the 
resultant determinant of the over-all system must van- 
ish and this gives the value of the reflection coefficient 
directly. Alternatively, A, and C, may be eliminated by 
means of (40) and (43), leaving a system of two sets of 
equations involving B, and D,. The required solution for 
any particular case is obtained in a straightforward 
manner, but the general details are too lengthy for 
inclusion here. 


measured values are also plotted in the above figures 
and in all cases are within the estimated experimental 
error from the computed values. The measured values 
were obtained by the usual tangent method, 7.e., by plot- 
ting the field minimum position in the partially filled 
guide vs short circuit position in the empty guide and 
subsequent analysis of the resultant curve. The dielec- 
tric slab was located in the slotted standing-wave de- 
tector section. 

It is interesting to note that the modulus of the re- 
flection coefficient is within one or two per cent of what 
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Fig. 3—(a) Modulus of reflection coefficient, (b) phase angle 
of reflection coefficient. 
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Fig. 4—Phase angle of transmitted wave. 


one would compute by assuming that the junction is 
equivalent to a junction of two transmission lines with 
characteristic impedances proportional to the respec- 
tive wavelengths in the empty and partially filled 
guides. From Fig. 3(b), it is seen that for values of 
t/b<0.82, the phase angle of the reflection coefficient is 
greater than mw radians, while for values of t/b>0.82 
the phase angle is less than 7 radians, corresponding re- 
spectively to more electric energy than magnetic energy 
stored in the evanescent modes at the junction and vice 
versa. The evanescent LSE modes store more magnetic 
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energy, while the evanescent LSM modes store more 
electric energy, except when the modes are close to 
propagating, so that 


2 2 


Tv T 
Yn? <— and f,?<—- 
a Ga 
For values of 
t 1 
= > 0.85, 1 res 
b a 


and the first evanescent LSM mode stores an excess of 
magnetic energy at the junction, resulting in a phase 
angle of less than 7 radians for the reflection coefficient. 

These results are obtained by considering only a 
small number of modes. If a larger number of modes are 
taken into account, sufficient compensation may take 
place so that the phase angle of the reflection coefficient 
does not become less than 7 radians for this particular 
sample. Further calculations are required to clarify this 
behavior. 

The modulus of the transmission coefficient can be 
computed only when the characteristic impedances of 
the empty and partially filled guide have been specified. 
Since these may be specified in any convenient way, the 
modulus of the transmission coefficient is not unique, the 
only restriction being that the transmitted power must 
be equal to the difference between the incident and re- 
flected power. 


CONCLUSION 


The use of the R-R method for obtaining approxima- 
tions to the first N eigenfunctions in a partially filled 
guide permits one to evaluate,in a straightforward man- 
ner, the junction discontinuity existing between an 
empty and partially filled guide. The reduction in com- 
putational labor is considerable, and the solution of 
transcendental equations and the evaluation of many 
complex expressions is avoided. The method outlined 
here may equally well be applied to the evaluation of 
the parameters of a slotted dielectric interface in free 
space. 


APPENDIX 


Eq. (32) in the text may be written in matrix form as 
follows 


ier Poi * > Ton Ba 
ee ‘ - Tyn ee 

- Pon doK 

=S YK Q : (44) 


anNK 
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or more briefly as 

[Tj ]laix] = vx?[Pii][eix] (45) 
and for the Ath solution as 

[Tis] [air] = ve? [Ps] [ae (46) 
Take the transpose of (2) to get 

{asx} [Tis] = vx? { esx} [Pic], 
Or 

Laue} [Tei] = vx? {esx} [Pas (47) 
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since 7;;=T;; and P;;=P;; and where fay} is a row 
matrix. Postmultiply (4) by [ajr|, premultiply (3) by 
fa;,.}, and subtract to get 


Cyn" — YR) (aix} [Pil lain | = 0. 


When yx?#/7r’, (5) gives 


N N 
yp S$ Onk@sRP sn = O, 
n=0 s=0 


(48) 


upon development of the matrix product. This proves 
the orthogonality of the eigenvectors. 


Coupling Through an Aperture Containing an 


Anisotropic Ferrite” 
DONALD C. STINSONT 


Summary—Coupling through an aperture containing anisotropic 
ferrites is investigated theoretically by asimple extension of Bethe’s 
small-hole coupling theory to include the dipole moment of the body 
in the aperture. The magnetic dipole moment of the ferrite body is 
ordinarily a vector but becomes a tensor upon the application of a 
magnetostati: field. This new theory is applicable to any situation 
where Bethe’s small-hole coupling theory is valid. Experimental 
verification was quite satisfactory and was obtained on two Bethe- 
hole type couplers: one with the waveguides parallel, and the other 
with the waveguides perpendicular. 


INTRODUCTION 


HE THEORY of coupling through small windows 

was formulated by Bethe more than a decade 

~ ago. Initially, he found that the amplitudes of the 

modes excited in a waveguide by a window were pro- 
portional to 


{zB xX He-rids 


where field 1 is the excited field, field 2 is a normal mode 
of the guide, and # is the inward normal. Later, he 
evaluated the integral over the window by developing a 
lumped-constant theory? for small windows and then 
applied this lumped-constant theory to side windows? 
in waveguides. 


* Manuscript received by the PGMTT, November 7, 1956. This 
work was supported by the U. S. Navy at the Univ. of Calif. under 
contract N7-ONR-29529 and is based on a thesis submitted in partial 
fulfillment of the requirements for the Ph.D. degree, Dept. of Elec. 
Eng., Univ. of Calif., 1956. 

+ Lockheed Aircraft Corp., Sunnyvale, Calif. 

1H. A. Bethe, “Formal Theory of Waveguides of Arbitrary Cross 
Section,” M.I.T. Rad. Lab. Rep. 43-26; March 16, 1943. 

2H. A. Bethe, “Lumped Constants for Small Irises,” M.I.T. Rad. 
Lab. Rep. 43-22; March 24, 1943. 

3H. A. Bethe, “Theory of Side Windows in Wave Guides,” M.I.T. 
Rad. Lab. Rep. 43-27; April 4, 1943. 


Bethe’s coupling theory depends upon his lumped- 
constant theory for small windows, which in turn de- 
pends upon replacing the excitation caused by the 
window by a quantity which is proportional to the fol- 
lowing parameters: 1) frequency; 2) the normal electric 
or tangential magnetic field (exciting field) which would 
exist at the center of gravity of the window if the win- 
dow were replaced by a solid metal wall; 3) the corre- 
sponding fields (induced fields) of the normal modes 
which are excited by the window; and 4) lumped con- 
stants (polarizabilities) which are functions only of the 
shape and dimensions of the window. The basis of his 
tumped-constant theory depends upon the fact that the 
excitation of the window can be replaced by “equiva- 
lent” electric and magnetic dipole moments. These 
“equivalent” electric and magnetic dipole moments lead 
him to consider the polarizabilities (which are defined as 
the “equivalent” dipole moments per unit incident 
field) as the true lumped constants of the window. This 
is logical since a window may act as either an inductive 
or capacitive element, depending upon its location and 
the propagating mode in the waveguide. 

Since his coupling theory applies only to cases where 
the window and the waveguides are filled with the same 
isotropic and homogeneous material, it is the purpose of 
this paper to extend his theory to include cases where 
the window is completely filled with an anisotropic fer- 
rite. The ferrite involved is anisotropic in the sense that 
its permeability becomes a tensor upon the application 
of a magnetostatic field. This extension will be made by 
adding the “equivalent” magnetic dipole moment of the 
ferrite to that of the window. 
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DriroLeE MOMENTS 


The ferrite in the window will affect both the electric 
and magnetic dipole moments. It is assumed that the 
isotropic permeability of the ferrite is real and unity 
for microwaves, and that the dielectric constant is dif- 
ferent from unity. The anisotropic permeability of the 
ferrite is a tensor and is given by expressions in the 
Appendix. Further, it is postulated that the dipole mo- 
ment of the filled aperture can be replaced by two other 
dipole moments: one to account for the aperture itself 
and the other to account for the material in the aper- 
ture. From Bethe?‘ it is known that the “equivalent” di- 
pole moments of the empty window are given as 


=i, = MM Hol ae MeHonm 
®., — Ey@n- Eo 


(Ja) 
(1b) 


where € is the electric inductive capacity of free space; 
Eo, Hy are the fields that would exist at the center of 
gravity of the window if the window were replaced by a 
metal wall; and 9%, Me, and @ are the magnetic and 
electric polarizabilities, respectively. 

Since the material in the window introduces a mag- 
netization and a polarization, assume that the “equiva- 
lent” dipole moments of the body in the window are 
given as 


Il, = NM ol a MWeM omm 
©, = On-P. 


(2a) 
(2b) 


Eqs. (2a) and (2b) are deduced by noting that, in 
general, the electric or magnetic dipole moment of a 
material is the product of its volume and the polariza- 
tion or magnetization. Since the polarizabilities 9%, 3%, 
and & have the dimensions of a volume and depend upon 
the shape of the window, it is assumed that the volume 
of the material can be replaced by the polarizabilities of 
the window when the material is located in the window. 

The total “equivalent” magnetic dipole moment of 
the filled aperture is the sum of (1a) and (2a); the total 
“equivalent” electric dipole moment of the filled aper- 
ture is the sum of (1b) and (2b). Thus, the following 
expressions are obtained for the total “equivalent” mag- 
netic and electric dipole moments: 


Mi Botl + M2Bontm 
Pn: Do 


— poll 
6 = 


(3a) 
(3b) 


where 

Bot = pol Lo Ty M,) 

Dot = 90 Eo. 
The quantity Q is a number greater than unity and ex- 
presses the fact that the polarization increases the mag- 
nitude of the electric dipole moment. The value of Q is 


determined from experiment and the quantity H)— Mo 
is determined from the theory of the anisotropic mag- 


4 Bethe, footnote 2, see (18) and (25). 
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netization of a ferrite. This is considered briefly in the 
Appendix. However, the result is that 


Boit == bitH oy =F Bint Hom 


Bont = pmitHor + bmmHom (4) 
where 
Kut = Lmnmt = bo(1 a X11) 
Bimt = — Bmit = — poXim. 


The magnetostatic field is applied in the m direction, 
while the microwave field is applied in the plane normal 
to n. 

GENERAL COUPLING EXPRESSIONS 


In this section, the expressions for the amplitudes of 
the normal modes coupled by a window between two 
waveguides is derived. Since this theory is well known, 
the presentation is brief. 

Silver® gives the following expressions for the field 
components of freely propagating modes in waveguides 
of arbitrary but uniform cross section: 

TE waves: 

Nab ae pale exp (+ 7802); ies 

Ey, = Eat exp (#79802); Bat 

Hi he» Br exp (+ 7802); Hat 
TM waves: 


K (wp) Wa 
ViWa x tz 
Balwu) 'Vaba. (5) 


I 


Ez,= + Gilsae exp (+ 7Ba2); Faz = Tie Oc 
LE; = Eat exp (+ 7B.z) ; EA = Vida 
Higie ct gs exp (+7842); Hat = weBa liz X Vida (6) 
Also, 
1 1 a aks 
Poo = = Sa= = f Ear X Mas tads 
D 2 
Pape) 


The functions Hy, Ea:, Eaz, Ha, are all real. Further, 
B, is the phase constant; w is the angular frequency; yu 
and € are the magnetic and electric inductive capacities, 
respectively; K is the eigenvalue of (V,?+K,?)F=0, 
where F is Wa or @a as the case may be; and w2ue= K,? 
+£6,”. Exp (jwt) time dependence is assumed and mks 
units are used throughout. The subscripts a and 6 in (5) 
and (6) designate the pair of mode indices mn and 
should not be confused with the transverse dimensions 
of rectangular waveguides. 

Also needed are expressions for the fields (5) and (6) 
when the waveguide is shorted at an arbitrary location: 


Short at g=d, —»<s<d: 
H, = 2H, exp (—jBad) sin Ba(z — a) 
E, = 2jEaz exp (—jBad) cos Balz — d) 
Ey = — 2j Ba: exp (—jBad) sin Ba(z — d) 
H, = 2Ha: exp (—jBad) cos Ba(z — d). (7) 


5S, Silver, “Microwave Antenna Theory and Design,” McGraw- 
Hill Book Co., Inc., New York, N. Y., art. 7.3; 1949. 
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Short at z= —d, —d<s<o:; 
H, = — 2Hqz exp (—j®ad) sin Ba(z + d) 
E. at aaa 24 Eas exp (—7B 4d) COs Balz + d) 
FE, = 2jEa: exp (—7Bad) sin Ba(z + d) 
Hi = — 2Hat exp (—jBad) cos Ba (z + d). (8) 


It has been shown by Bethe! and by Silver® that the 
modal amplitudes of the fields set up in an infinite 
waveguide by a window in the wall are the following: 


2A wa 


= if Ty > H.-- nds 


WwW 


(9) 


2B eoD in 


= f E, X Het- ads. (10) 


WwW 


The + and — superscripts on the field H2 indicate waves 
going in the direction of positive or negative z, respec- 
tively. The field F, is the electric field set up in the 
guide by the window. The field Hz is a normal mode 
field of the waveguide when the window is absent. The 
surface integral is over the window, and the set of axes 
l, m, n are fixed in the aperture. Further, the / direction 
is always parallel to the long dimension of the aperture, 
the m direction is always parallel to its narrow dimen- 
sion, and m points away from the source of excitation. 

Eqs. (9) and (10) are also valid when the excited 
waveguide is semi-infinite and when either side-window 
or end-window coupling (iris coupling) is being con- 
sidered. In either case, it is only necessary to obtain the 
field Hz in (9) and (10) from (7) or (8), rather than from 
(5) or (6). 

The integrals in (9) and (10) have been evaluated by 
Bethe? when the window is empty. When the window is 
filled with a ferrite, they become 


2A aSa . Leet 
= — jo( FIM BoitHe + M2Bom*Hom 


(11) 
I) Sy ON 
ak CDont Eon). 


The fields Bot and Dot are defined in (3) and (4). The 
upper signs in (11) refer to A,; the lower signs refer to 
B,. For a circular window of radius 7, 1% = M2 = (4/3) 
and @=(2/3)r*. For apertures of other shapes, S% al- 
ways corresponds to an incident magnetic field parallel 
to the long dimension of the aperture, while 2 refers to 
the narrow dimension of the aperture. It should be noted 
that 9% and MM. are functions only of the shape of the 
window and are not to be confused with the fields 1 and 
2. For a specific application, the fields E., Hy will be 
defined by one of the expressions (5)—(8); the fields Eo, 
Hy will be defined by one of the same expressions except 
that the mode index a is replaced by 8. Eq. (11) is also 
valid when the waveguides are semi-infinite and for iris 
coupling between waveguides. 


6 Jbid., art. 9.10. : 
7 Bethe, footnote 2, see (51) and (55). 
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CouPLING THROUGH SIDE WINDOWS 


In this section, the general expression (11) is evalu- 
ated for two particular cases: 1) the axis of the primary 
guide (source of excitation) is parallel to the axis of the 
secondary guide (excited guide), and 2) the axis of the 
primary guide is perpendicular to that of the secondary 
guide. The first device is designated as a parallel 
coupler,? and the second device as a perpendicular 
coupler. The set of axes fixed in the primary waveguide 
is denoted by &, 7, €; the set of axes fixed in the second- 
ary waveguide is denoted by x, y, z. As mentioned be- 
fore the set of axes /, m, n are fixed in the aperture. The 
orientation of the three sets of axes is illustrated in 
Fig. 1. However, m and z are made parallel for the two 


Fig. 1—Orientation of axes. 


couplers under consideration. Also, the mode index a 
applies to the secondary guide, and the mode index 6 
applies to the primary guide. Case A is defined as the 
situation where identical rectangular waveguides are 
joined on their broad sides, the TEi») mode is propagated, 
and the window is centered. 


8 Both of these waveguide configurations were considered by A. D. 
Berk and E. Strumwasser, “Ferrite directional couplers,” Proc. IRE, 
vol. 44, pp. 1439-1446; October, 1956. However, their work was based 
on the theory of scattering by an obstacle in a waveguide. Further, 
they considered only ferrite cylinders extending into both wave- 
guides with the coupling holes located at positions of circular polar- 
ization of the magnetic field. The theory presented in this paper is 
also applicable to these situations and is sufficiently different, it is 
felt, to warrant a paper on these configurations. Such a paper is 
nearing completion. 

® This particular waveguide configuration was also studied by 
R. W. Damon, “Magnetically controlled microwave directional 
coupler,” J. Appl. Phys., vol. 26, pp. 1281-1283; October, 1955. His 
work is of a qualitative nature and considers a ferrite cylinder located 
at the position of circular polarization of the magnetic field. Further- 
more, his theory is also based upon an extension of Bethe’s coupling 
theory. However, his theory is not of a general nature but merely re- 
places the magnetic dipole moment of the hole by the magnetic dipole 
moment of the ferrite body in the hole. 

It would seem that we conceived the idea of extending Bethe’s 
coupling theory about the same time. The author originally felt that 
Bethe’s method could be extended to treat the case of an anisotropic 
ferrite filling the coupling hole (Inst. of Eng. Res., Progress Rep. 
Ser. No. 60, Issue No. 7, Electronics Res. Lab., Univ. of Calif., 
Berkeley, Calif., p. 11; January 15, 1955) and later developed such a 
theory during the summer of 1955 (Progress Rep. Ser. No. 60, Issue 
No. 10, p. 9; October 15, 1955). 
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Parallel Coupler, Primary and Secondary Waveguides If Bof from (15) and Do from (3) are substituted into 
Infinite (11) after using the proper field components from (1) or 


Expressions for A, and B, are obtained from (11), (2), the following expression results: 


where the fields £2, H, are specified in (5) or (6). The 24,5, 

fields Dot, Bot are given in (3) and (4); the fields Eo, Ho 

are given by (1) or (2). This results in 

DEA Sa 

eBada 

—— jo |(F GA ant + JM eH asin) Aor + 7( FMA a ipimt 
a INGE aphinnt iL + €o POE anE tn |. (12) 

It should be noted that the transverse components of 


the fields with the mode indices a and 6b have been left 
in the aperture coordinates, since all that was specified 2BaSa 


thus far _ that mand z and ¢ are parallel. However, once For a round coupling hole of diameter d, (17) simplifies 
the orientations of the /m, én, and xy axes are prescribed, 


= 70) Li( +9 oH amrit —M eH 4¢H azbmnt) 
DB Sy 
et M1 A oA aint + MWe o¢H azmit + €9 POEs, Ean |. (16) 


With the field components in (16) given by (5), the 
following expression results for Case A: 


24 oSa 
= — jo[ FM11?6?(«o2y2a2)— wnt + eom?a~2PQ |. (17) 


. to 
one should insert the field components in their proper 
coordinate systems. Therefore, if (5) is used in (12), the Aq 
following expression results for Case A: = —jrd*(3abd 9) Fal Fx yt1/2(d¢/do)?OF af? ]. (18) 
DA Sa Ba 
ee & cui Perpendicular Coupler, Primary Waveguide Semi-Infi- 
=.— NG 2R2 22s 1 if 2 2 : ; i { 
jo| F Man's (otura?) wut + eom?a 0]. (13) nite, Secondary Waveguide Infinite 
2 Beda 


This situation is the same as the preceding case, ex- 
In (13), use has been made of the fact that Ha,= —H,,, cept that the fields H;, H;, and E, are given by (7) with 


ep pe) am = ag ANG im = ey: the mode index a replaced by 8, instead of by (5) or (6). 
If a round coupling hole of diameter d and thickness ¢ Thus, for a short at z=d, — «© <z<d, and window at 
is considered, (13) reduces to the following: the origin: 
AeSe 
= — jw exp (—jBrd) | —( FMA aipimt + JM2H actimmt) Hoe cos Bod (19) 
Bie 
= (+9. Hamat on IM eH aclimt®) A o¢ ie jeoPQEanE sn | sin Bd} : 
Ag Using the field values from (5) in (19), the following re- 
= —jrd3(3abd,) Fx [ F (1 — x22) +1/2(Xg/do0)?OF af x *|(14) lation is obtained for Case A and 6,d=q7, where g is an 
: integer: 
’ AgSa = — BaSq = jwoMir?B?(w7p72a?) wit. (20) 


where Fy=Fa=exp{ —2mtd.-1[1—(Ac/do)?]/2}; Ac is 
different for Fg and Fy,!° since the electric and mag- For around coupling hole of diameter d, (20) reduces to 
netic fields (although below cutoff) in the window propa- : 

; A= —. By, = — 427d* (Babys) FP xxoy 2 
gate as different modes. It was assumed that the terms L2EE CIN iene! (21) 
Fz and Fy are not affected by the presence of the ma- ExprrIMENTAL RESULTS FOR SIDE-WINDOW COUPLING 


ene aang aes In the last section, the amplitudes of the modes 


Perpendicular Coupler, Primary and Secondary coupled by the parallel and perpendicular coupler were 
evaluated. Here, the coupled power for Case A condi- 
tions and for a round coupling hole is calculated. The 
power coupled into the secondary waveguide by the 
window is equal to the square of A, or Ba, since these 
amplitudes have been derived for unit amplitude inci- 
dent fields in the primary waveguide. 

Bot = — wits — bint, The coupled power for (14) becomes 


Waveguides Infinite 

Expressions for A, and B, are again obtained from 
(11), and the fields E2, Hz are specified in (5) or (6). 
However, the fields Bo# must be made proportional to 
the proper fields in the primary waveguide. Thus, 


Bont = — pms — pmb Hy. (15) cyt = Cy + 10 log { [+(1 — C'A) 


10 C, G. Montgomery, “Technique of Microwave Measurements, ” 


McGraw-Hill Book Co., Inc., New York, N. Y., p. 862; 1947. + 1/2(X4/do)2?OF Fu]? + (C’B)?} (22) 
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where 
Gy = 20 log [rd2(3abd 4) Wy |. 


The upper and lower signs refer to the mode coupled in 
the same and opposite directions, respectively, as that 
of the incident mode. The quantities A, B, and C’ are 
defined in the Appendix. For d=0.1235", 7020204 
Fy=0.555, and f=9.350 kmc, the unperturbed mag- 
netic coupling (Co) is —55.3 db. Thus, (22) becomes 


— 55.3 + log {[#(1 — C’A) + 0.81890]? 
+ (C’B)*}. (23) 


Cyt = 


When no magnetostatic field is applied, the ferrite is 
isotropic, and C’ vanishes. Consequently, the factor Q 
can be evaluated experimentally from the zero mag- 
netostatic field expressions. The value used here in com- 
paring theory and experiment will be an average value 
rather than one which gives the best agreement. There- 
fore, the resultant curves will not agree exactly at the 
zero magnetostatic field point, but this is acceptable 
since the interest here is in qualitative agreement. 

Comparison of the theoretical expression (23) with 
experimental results is offered in Fig. 2 for a Ferramic 
A sphere. (Parameter Mi, is proportional to the damping 
constant and is defined in the Appendix.) According to 
theory, the diagonal susceptibility is an even function of 
the applied magnetostatic field. Since the experimental 
curves for both directions of magnetostatic field were 
the same within the limits of experimental error, only 
the experimental curves for one direction of the applied 
magnetostatic field is presented. 


Magnetostatic field (Kilo oersteds) 


Fig. 2—Forward and reverse coupling in parallel coupler, 
Ferramic A sphere, both waveguides infinite. 


The effect on the reverse coupling of a variation in Q 
is shown in Fig. 3 for a Ferramic A sphere. Note that an 
increase in Q increases the initial coupling slightly and 
also reduces the variation between the maximum values 
of coupling. This is correct since an increase in Q corre- 
sponds to an increase in the dielectric constant, which 
in turn corresponds to an increase in electric coupling. 
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Magnetostatic field (kilo oersteds) 


= 740) 


Fig. 3—Reverse coupling in parallel coupler as function of dielectric 
properties of Ferramic A sphere, both waveguides infinite. 


Magnetostatic field (Kilo oersteds) 


Fig. 4—Reverse coupling in parallel coupler as function of saturation 
magnetization of Ferramic A sphere, both waveguides infinite. 


Fig. 4 shows the effect on the reverse coupling of a 
variation in saturation magnetization. A reduction of 
this parameter causes a reduction of the coupled mag- 
netic field, thereby reducing the perturbing effect of the 
ferrite. 

The third parameter of the ferrite studied is the 
damping constant. The behavior of the reverse coupling 
for various reduced damping constants is shown in 
Fig. 5. It should be remembered that the reduced damp- 
ing constant is the ratio of the actual damping constant 
to the magnetization (saturation magnetization here). 
This parameter is the most critical of the three, as one 
would expect in aresonance-type phenomenon. Note that 
a reduction of the damping constant causes an increase 
of the perturbing effect of the ferrite. 

The coupled power for (18) becomes 


— 49.4 + 10 log [(+ EF + 0.8450)? 
+ (BG)?] (24) 


where d=0.1495 inch, ¢=0.020 inch, Fy-0.617, f =9.350 
kmce, and the unperturbed coupling (Cy) is —49.4 db. 
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The quantities EF and KG are defined in the Appendix 
and are odd functions of the applied magnetostatic field. 
This means that a reversal of the direction of the mag- 
netostatic field results in an interchange of the forward 
and reverse couplings: that is, C+(+)=C-(—) and 


Magnetostatic field (Kilo oersteds) 


(e) | 2 3 4 5 
fe) 
Q= 1.33 
C'= 9045 
-10 = a 
-20 
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Fig. 5—Reverse coupling in parallel coupler as function of reduced 
damping constant of Ferramic A sphere, both waveguides infinite. 


C+(—)=C-(+), where, for instance, C+(—) indicates 
the forward coupling for negative values of the applied 
magnetostatic field. Comparison of theory and experi- 
ment for positive values of the applied magnetostatic 
field is illustrated in Fig. 6 for a Ferramic A sphere. The 
agreement between theory and experiment is accepta- 
ble, although it is felt that better quantitative agree- 
ment can be obtained by choosing a smaller value for the 
reduced damping constant. 
The coupled power for (21) becomes 


Cit = — 43.4 + 10 log [£°(F? + G’)] (25) 
where Cy is 6 db larger because of the short in the 
primary waveguide. The constants of the window are 
the same as for (24). Since the short in the primary 
waveguide is located so as to annul the electric coupling, 
there is no coupling into the secondary waveguide until 
the ferrite becomes anisotropic. Note in Fig. 7, for Fer- 
roxcube 4-A, that a small value of the magnetostatic 
field is sufficient to increase the coupling to the value 
it would normally be when only the electric field couples. 
Agreement is quite satisfactory except for small values 
of applied magnetostatic field. One sees very clearly 
that the theoretical curve predicts too large a coupling 
initially. This is probably caused by the fact that the 
actual magnetization has been replaced by the satura- 
tion magnetization. The agreement at resonance could 
also be improved by using a smaller reduced damping 
constant. 

Although the theory derived is valid only when the 
coupling aperture is completely filled with a material, 
it is of interest to examine experimentally the effect on 
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Fig. 6—Forward and reverse coupling in perpendicular coupler, 
Ferramic A sphere, both waveguides infinite. 


Magnetostatic field (Kilo oersteds ) 


Db 


= 109) 
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Fig. 7—Reverse coupling in perpendicular coupler as function of re- 
duced camping constant of Ferroxcube 4-A sphere, primary wave- 
guide semi-infinite, secondary waveguide infinite. 


the magnetic coupling of partially filling the aperture 
with a ferrite. In this case, shorts were used to annul 
the electric coupling. However, the electric coupling could 
also be excluded by properly choosing the location of 
the window. For comparison purposes, the magnetic 
coupling curves for a spherical, a rectangular, and a 
disk-shaped Ferroxcube 4-A sample are shown in Fig. 8. 
The spherical and disk-shaped samples were placed in 
a round aperture; the rectangular sample was placed 
in a rectangular aperture with the long dimension 
parallel to the axis of the primary waveguide. These 
were all obtained with the perpendicular coupler and for 
Case A conditions. The curves for the rectangular and 
disk-shaped samples are very similar. However, the 
curve for the spherical sample has a much larger ampli- 
tude at resonance than either of the other curves, and 
its resonance also occurs for a much smaller value of 
applied magnetostatic field. 
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Fig. 8—Effect on magnetic coupling in perpendicular coupler by 
partially filling aperture with various shaped Ferroxcube 4-A sam- 
ples. Jf 0.0965 inch dia. sphere in 0.151 inch dia., 0.020 inch thick 
aperture. A 0.250 inch X0.075 inch X0.020 inch thick sample in 
aperture of same size and 0.010 inch thick. © 0.100 inch dia., 
0.021 inch thick disk in 0.151 inch dia., 0.020 inch thick aperture. 


CONCLUSION 


The theory presented is satisfactory for small samples, 
although the sample sizes used here are near the upper 
size limit for these materials. For materials with smaller 
losses, it would be necessary to use even smaller sample 
sizes in order to avoid dimensional effects. 

The theory of the magnetization of the ferrite is the 
weak link in this coupling theory. Consequently, any 
effort to apply this coupling theory to larger samples 
would require a proper modification of the magnetiza- 
tion expressions. It is felt that the dimensional effects 
in the sample enter in two ways. The first occurs because 
of the inhomogeneity in magnetization and can be ac- 
counted for by properly choosing the constants of the 
ferrite. In other words, we can choose values for the con- 
stants in such a manner as to obtain good agreement be- 
tween theory and experiment but these values will not 
necessarily be the true constants of the material. The 
other effect occurs when the sample becomes electrically 
large and acts as a resonator. The present coupling the- 
ory does not account for this. 


APPENDIX 
TENSOR MAGNETIZATION OF A FERRITE 


The purpose here is to evaluate the susceptibilities 
defined in (4). These susceptibilities result when a fer- 
rite is made anisotropic by the application of a magneto- 
static field in the direction. Although this has been 
done by several authors,'!!” part of the work is repeated 
here for continuity. The form followed is that given by 
Beljers.!! 


11H, G. Beljers, “Measurements on gyromagnetic resonance of a 
ferrite using cavity resonators,” Physica, vol. 14, pp. 629-641; Feb- 
rurary, 1949. 

12°C, L. Hogan, “The microwave gyrator,” Bell Sys. Tech. J., vol. 
31, pp. 1-31; January, 1952. 
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The fundamental equation for the magnetization, M, 
is the following: 
My =MxB-\B-BMM-M) (26) 
where y=—ge/2m (the magnetomechanical ratio, a 
negative quantity for an electron); e and m are the elec- 
tron charge and mass, respectively; g is the spectro- 
scopic splitting factor; \ is a damping constant; and 
B=u(H+M). The damping constant can be intro- 
duced in different manners, but the two commonly used 
forms both give identical results for the order of ap- 
proximation used here." 
For the problem under consideration, the following is 
chosen: 


H, hy — (Ni — 1/3) Mi 
lal. po ie. a Nes — 1/3)Mm 
He = He(Nnaet tM 


I 


where the N;’s are the demagetizing factors, the quan- 
tities h; and h, are the applied microwave magnetic 
fields, and H is the applied magnetostatic field. 
Consistent with the usual assumption, the microwave 

magnetizations are as follows: 
My, = yML[hi{ wy + 2) [HZ + M(Nin — Nz) ] 

a jopuodr} = jopolm | 
Mm = yML-[hm{ wory(1 + 2) [H + M(Ni — N,)] 

— jopods} + joouohi] (27) 


where 


L= = wo + po’y2(1 + 2) [HH + MM. — Nn 
-[H + M(Nn — Nm) | — 2jopoyH 
— jopoyaM (Ni — Nm — 2Nn), 
\i1=A/WM is a reduced damping constant, and M is the 
saturation magnetization in the m direction. When the 


ferrite sample is spherical, the magnetizations simplify 
considerably to the following: 


Mi = xuhi t+ xXimhn 


My, = Xmily He Nneaias (28) 


where 
XI eX ta IX ui 
Xim = Ximr oh Ihas 


For computational purposes, the following substitutions 
are made: 
rey 


BG boy Hw (1 =F Na2) 22, (29) 


__'8 We introduce our damping constant according to Landau and 
Lifschitz (see Beljers, op. cit.). 
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Susceptibility 


Xgq = CCA + JB) 


Ch = Yt Me oMw~! C1 +A2)4 


t 


x= Wotw “Cl +d 2) 


ae ae 


— | = 0.05 


) 0.4 0.8 1.2 1.6 2.0 
Normalized magnetostatic field (Xx) 


Fig. 9—Diagonal susceptibility of small ferrite sphere as 
function of normalized applied magnetostatic field. 


Thus, the components of the susceptibilities become 


Xir = C’A 
Xui = C'B 
Ximr = EF 
Ximi = EG (30) 


where 


Cl = x’ uoMo "(1 + dy?)1/? 


AER AD) 

B= —a'/b 
E = y'wMoa-! 
Fi= cf 
G=-—¢/b 


/ 


x { a? — [1 — 2r12(1 + 2)-1)} 
a’ = dy(a? + 1)(1 + AW2)71? 
6 de x 
6 = Wield + AY)? 
b = x4 — 2?[1 — 221 + A)1] 4 1. 
Note also that Xmm=x1u, and Xim = —Xmi- 


In most ferrite work the susceptibilities are defined so 
that 


Mi 
Mn 


xhi — 7K hn 


jKhi ain Xlom (31) 
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Xpm = ECF + JG) 


Susceptibility 


(@) 4 8 ka 1.6 2.0 
Normalized applied magnetostatic field (x) 
Fig 10—Nondiagonal susceptibility of small ferrite sphere as 
function of normalized applied magnetostatic field. 


where 


Y= x’ — jx” 


K= Ko GK", 
Using this notation, 
x =C'A 
x = —C'B 
K’ = — EG 
VG Sa tsi (32) 


For the simple case of no damping, the imaginary parts 
of both x and K vanish. Eq. (33) then reduces to the 
one that is characteristic of a gyrotropic medium. 

In order to understand better the behavior of the 
permeability of the ferrite when it is anisotropic, several 
families of curves of x1: and xim are plotted as a function 
of x, with M and ); as parameters; Fig. 9 is a plot of 
xu vs x for \,=0.05 and \,=0.10. Similarly, Fig. 10 
is a plot of Xm vs x for \:=0.05 and \,; =0.10. In both 
figures, the coefficients C’ and E were set equal to unity. 
It should also be noted that x1 is an even function of the 
applied magnetostatic field, whereas xim is an odd func- 
tion of the applied magnetostatic field. 
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An Adjustable Sliding Termination for 
Rectangular Waveguide" 


ROBERT W. BEATTYT 


Summary—A new adjustable sliding termination for rectangular 
waveguide has been developed. The termination is of simple design 
and can easily be adjusted to have reflection coefficients from zero 
to nearly unity in magnitude and any desired phase. In addition to the 
usual applications of adjustable sliding terminations for rectangular 
waveguide, it provides a suitable design for an adjustable transfer or 
secondary standard of impedance for rectangular waveguide systems. 


INTRODUCTION 


DJUSTABLE sliding terminations for rectangular 
waveguide are useful for a variety of purposes, 
including calibration of microwave impedance 

measuring devices, obtaining the directivity of direc- 
tional couplers, and determining the parameters of 
waveguide junctions. 

An adjustable sliding termination for rectangular 
waveguide was described by Grantham in 1951.! Sev- 
eral years later, a similar termination was produced 
commercially. This termination? differed from its prede- 
cessor in details of construction of the dissipative ele- 
ment and the reflecting antenna. Both of these termina- 
tions were designed primarily to provide minimum re- 
flection and could not be adjusted over a wide range of 
vswr (voltage standing-wave ratio). The commercial 
termination has a vswr range of 1.005 to 1.15. The 
principle of the double slug tuner was used* to obtain 
an adjustable sliding termination with a somewhat 
greater range of adjustment. An adjustable sliding 
termination having a wide range of vswr was described 
by Kato and Sakai,‘ but it was not possible to adjust 
this termination for cancellation of reflections. 


PRINCIPLE OF OPERATION 


As shown in the diagram of Fig. 1, this termination 
slides inside a rectangular waveguide and consists of a 
short-circuiting piston to which is attached a dissipative 
strip supported by a dielectric rod, which can rotate and 
slide relative to the piston. The phase of the reflection 
from the strip can be varied by sliding the strip, while 
the magnitude of the net reflection from the short cir- 


* Manuscript received by the PGMTT, November 19, 1956. 

{+ Natl. Bur. of Standards, Boulder, Colo. 

1R. E. Grantham, “A reflectionless wave-guide termination,” 
Rev. Sci. Instr., vol. 22, pp. 828-834; November, 1951. 

2W. A. Andrews, U.S. Patent No. 2,701,861. 

3R. C. Ellenwood and W. E. Ryan, “A uhf and microwave 
matching termination,” Proc. IRE, vol. 41, pp. 104-107; January, 
1953. 

4N. Kato and T. Sakai, “Waveguide Type Variable Impedance 
Circuit and its Application for Rieke Diagram,” Rep. of the Micro- 
wave Communication Res. Committee in Japan, p. 7; December, 
1955. 


cuit can be varied by rotating® the strip. With inde- 
pendent control of these two motions, complete cancel- 
lation of reflections can be obtained. On the other hand, 
with the strip surface perpendicular to the electric field, 
minimum losses occur in the strip and almost perfect 
reflection is obtained. It is possible to adjust the termi- 
nation to any intermediate condition, then to slide the 
entire assembly, to obtain almost any reflection coef- 
ficient desirable. 


RECTANGULAR WAVEGUIDE 
DISSIPATIVE STRIP 


> 
RSSGHD 


a 
ORL 


DIELECTRIC ROD SHORT-CIRCUITING PISTON 


Fig. 1—Essential features of termination. 


THEORY 


Apparently, the analysis of the fields inside a rec- 
tangular waveguide containing an arbitrarily positioned 
dissipative strip has not yet been published. Lacking 
this basic information, a rigorous analysis of the action 
of this termination has not been attempted. An ap- 
proximate study based upon microwave circuit theory is 
given below. 

As shown in Fig. 2, the termination can be regarded 
as an attenuator terminated in a short-circuited line of 
variable length. In terms of the scattering coefficients® 
Su, Siz, and Sx» of the attenuator, the input voltage re- 
flection coefficient is 


S12” 


T; = Su as ata ea 
Seo + 6776! 


(1) 


where 8=27/g@, and Xg equals the wavelength in the 
waveguide. The condition for cancellation of reflections 
(T; =0) is 


S32? = Sy1e7'(1 + Soe 7784), (2) 


5 G. C. Southworth, “Principles and Applications of Waveguide 
Transmission,” D. Van Nostrand Co., Inc., New York, N. Y., pp. 
374-376; 1950. An attenuator employing a rotating strip is described. 

6 Scattering coefficients are the elements of the scattering matrix. 
See C. G. Montgomery, R. H. Dicke, and E. M. Purcell, “Principles 
of Microwave Circuits,” M.I.T. Rad. Lab. Ser., McGraw-Hill Book 
Co., Inc., New York, N. Y., vol. 8, pp. 146-149; 1948. 
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Fig. 2—Approximate equivalent circuit representation. 


Normally, [Su | and LSes | are much less than unity, and 
the approximate condition for cancellation of reflections 
is 


S19 Se? e% 


(3) 


Rotation of the strip will, in general, change both 
|Sie?| and |S |. A typical variation of these quantities 
with 0 (the angle between the electric field direction and 
the normal to the surface of the strip) is shown in 
Fig. 3(a). At the angle 6c, the magnitudes of Sj)? and Sy, 
are equal, and it is possible to obtain cancellation of re- 
flections by sliding the strip, varying /. If the strip is 
too short, |S1.2| will not decrease enough to equal [Si |, 
as shown in Fig. 3(b). However, it is possible to increase 
[Six | by adding a reflecting object at the end of the 
strip, so that cancellation of reflections can again be 
obtained as in Fig. 3(c). 


UNIT CIRCLE 8) 


fe 12 | \ c 2 
[Se— | | \ [Be 


(a) (b) 


Fig. 3—Dependence of |.Su| and |Si.?| upon rotation of thin rec- 
tangular strip: (a) rectangular strip 0.35 inch X 1.5 inches, (b) rec- 
tangular strip 0.35 inch X0.75 inch, and (c) metal reflector added 
to increase | Si]. 


DESIGN 


A number of considerations can influence the design 
of the termination. If the primary need is for a reflec- 
tion-free termination, with no need for a wide range of 
vswr, it seems advisable to use a thicker strip of dissipa- 
tive material such as Synthane. The extra thickness can 
give added strength to a long, tapered strip, which will 
require less rotation to achieve cancellation of reflection, 
and give a smoother adjustment which is less frequency 
sensitive than with shorter strips. Examples of such a 
design are shown in Fig. 4, and a diagram of typical 
variation of | Si2?| and | Sia is shown in Fig. 5. 

If a wide range of adjustment of vswr is desired, a 
thin strip is required. IRC resistance strip may be used 
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2 $2 3 : : 4 


Fig. 4—Adjustable sliding loads for different waveguide sizes. 


Fig. 5—Variation of |.Su| and |.Si.2] with rotation of 
thick tapered strip. 


to obtain fairly high ranges, but a dissipative film on a 
thin mica strip is capable of greater range. A thin film 
has little strength and it may be desirable to add a re- 
flecting disk or rod as shown in Fig. 6. The disk or rod 
permits a shorter strip to be used, while retaining the 
ability to cancel reflections. 

Control of the mechanical motion required may be 
achieved by the arrangements shown in Figs. 7 and 8, 
next page. Independent control of the sliding and rota- 
tion of the strip is achieved in both arrangements, but in 
Fig. 8, the strip may be rotated and slid by hand until 
the clamping screw is used, which clamps the dielectric 
rod supporting the strip so that it can only be moved by 
the limited fine adjustments A and B. One adjustment 
of sliding (knob A) and an independent adjustment of 
rotation (knob B) are provided. 

As a matter of practical interest, an adjustable load 
having satisfactory performance over the recommended 
frequency range of RG-52/U waveguide was con- 
structed, using a rectangular IRC resistance strip 200 
ohms per square, 0.35 inch X1.50 inches. The strip was 
not tapered, but mounted on a dielectric rod as was 
shown in Fig. 1. 
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Fig. 6—Short strip and reflecting disk. 


Fig. 7—Mechanical controls permitting independent adjustment 
of rotation and sliding, with control locking. 
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Fig. 8—Mechanical controls permitting coarse manual adjustment 
and locking, followed by fine independent adjustments of rotation 
and sliding. 
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Field Displacement Isolators at 4, 6, 11, and 24 KMC* 


S. WEISBAUMY anp H. BOYET{ 


Summary—Performance of ferrite field displacement isolators at 
various frequency bands is described. Single and double-slab isola- 
tors have been constructed in rectangular waveguide. Four single- 
slib isolators are reported in the following frequency bands: 3700-— 
4200 mc; 5925-6425 mc; 10,700-11,700 mc and 23,500-24,470 mc; 
one double-slab isolator is described in the frequency range 10,700- 
11,700 me. 


INTRODUCTION 


7 HE field displacement isolator consists of a rec- 
tangular slab of ferrite partially filling a rectangu- 
lar waveguide.'! The ferrite slab is transversely 

magnetized and has resistance material appropriately 

disposed on one of its sidefaces (see Fig. 1, opposite). 

Basically, this device produces isolation for the follow- 

ing reasons: When a microwave travels through the de- 

* Manuscript received by the PGMTT, December 5, 1956. 
+ RCA, New York, N. Y. Formerly with Bell Telephone Labs., 

Murray Hill, N. J. : 

t Bell Telephone Labs., Murray Hill, N. J. 
1 A. G. Fox, S. E. Miller, and M. T. Weiss, “Behavior and appli- 


cations of ferrites in the microwave region,” Bell Sys. Tech. J., vol. 
34, pp. 5-103; January, 1955. 


vice in a given direction, the polarization of the rf mag- 
netic field in the plane perpendicular to the biasing 
magnetic field is opposite in sense to that which it would 
be for the reverse direction of propagation. As a result, 
the rf magnetic field interacts differently with the spin 
precessions of the electronic magnetic moments for the 
two directions of propagation, and this leads to different 
electric field distributions across the waveguide for the 
two directions of propagation? (see Fig. 2). If this dif- 
ference in electric field strengths at one face of the ferrite 
can be made large, resistance material, suitably placed 
on the ferrite, will attenuate reverse traveling waves toa 
much larger degree than forward traveling waves, and 
isolation is thereby effected. 

The above principles have been applied with success 
to an isolator at 6 kmc (0.2-db forward loss, 30-db re- 
verse loss, 1.06 vswr over the band 5925-6425 mc). The 


_ ?B. Lax, H. J. Button, and L. M. Roth, “Ferrite phase shifters 
in rectangular waveguide,” J. Appl. Phys., vol. 25, pp. 1413-1419; 
November, 1954, 
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Fig. 1—Assembly of ferrite, resistance material, magnet, and wave- 
guide in single-slab ferrite field displacement isolator and disposi- 
tion of resistance material on ferrite slab. 
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Fig. 2—Theoretical transverse electric field distribution in 
single-slab isolator, full height ferrite. 
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experimental, theoretical, and design details for this 
device have been elaborated upon elsewhere.’ Subse- 
quent work with the field displacement isolator has met 
with success at other frequency bands: 3700-4200 mc; 
10,700—-11,700 mc; and 23,500—24,470 mc. The purpose 
of this paper is to summarize the performance of field 
displacement isolators developed to date in the various 
frequency bands at Bell Telephone Laboratories. 
Before doing so, it is desirable to discuss briefly two 
important features of the field displacement isolator: 
the partial height ferrite geometry used, and the dis- 
position of the resistance material. Less than full height 
ferrite slabs are used to insure good match into the 
isolator. With such partial height slabs the boundary 
requirements on the top and bottom faces of the ferrite 
are less stringent and a smaller portion of the incident 
wave needs to be converted into a reflected wave to 


3S. Weisbaum and H. Seidel, “The field displacement isolator,” 
Bell Sys. Tech. J., vol. 35, pp. 877-898; July, 1956. 
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satisfy the boundary conditions. This results in con- 
siderably smaller vswr for the partial height ferrite than 
for the full height ferrite. 

A disposition of the resistance material was conceived 
which makes use of the nonreciprocal nature of the 
longitudinal components of electric field set up by the 
partial height ferrite.’ Resistance material is placed in a 
region where both the fundamental transverse com- 
ponent and longitudinal component of electric field are 
small for the forward direction but, as a result of the 
nonreciprocity of both components, both may be large 
at the resistance material for the backward direction of 
propagation. This configuration of resistance material is 
shown in Fig. 1. The actual resistance material is a 
baked resistance coating of graphite and resin covered 
by a baked clear coating, similar to that used in attenu- 
ator vanes. The length of the resistance strips is critical 
in determining the forward and reverse attenuations. 
Experiments indicate that the attenuation is not a linear 
function of length of resistance strip. 

As an aid in the design of isolators at various fre- 
quency bands, we have relied to a first approximation on 
a scaling*® of transverse geometric quantities and mag- 
netic quantities. Thus starting from a frequency w, we 
may scale to a new frequency fw (f any number) by 
changing all waveguide and ferrite geometries and spac- 
ings by a factor 1/f and all magnetic quantities by a 
factor f (see Appendix). 


PERFORMANCE 


In Figs. 3, 4, 5, and 6, next page, we show graphically 
the performance of the field displacement isolator at 4 
kmce, 6 kmc, 11 kmc, and 24 kmce, respectively, together 
with the values of the operating parameters. 47M, de- 
notes the saturation magnetization of the ferrite. The 
other symbols are defined in Fig. 1. The length of ferrite 
is 5 inches in each isolator, but shorter lengths are cur- 
rently being investigated. 

In the 24-kmce case the largest 47M, ferrite available 
was 4900 Gauss. 

A double-slab field displacement isolator was con- 
structed at 11 kmc. In this construction the two slabs 
are located symmetrically with respect to the center line 
of the rectangular waveguide. The two slabs are mag- 
netized by equal but oppositely directed magnetic fields 
so that the spin precessions in each interact equally with 
the rf magnetic fields at each slab. Each slab has a re- 
sistance strip configuration identical with that shown in 
Fig. 1. The length of ferrite slabs is 5 inches. Perform- 
ance and operating conditions are shown in Fig. 7. 

Table I, p. 197, summarizes performance of most re- 
cent isolator models. All isolators are single slab except 
where noted. 


4S. Weisbaum and H. Boyet, “A double slab ferrite field displace- 
He isolator at 11 kmc,” Proc. IRE, vol. 44, pp. 554-555; April, 
1956. 
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ey ci eaten nee The theoretical basis of scaling to other frequencies is 
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presented in the Appendix. Table II opposite presents 
the geometrical and magnetic operating conditions to be 
expected in the various frequency band isolators, using 
avers cee the 6-kmc isolator as a basis and scaling the parameters 
Peps Beaabe up (or down) to the other frequency band _ isolators 
(4 kmc, 11 kmc, 24 kmc). The values in parentheses are 
the values expected from scaling, the others are the 
values used in the actual device. 
In all cases, the cross sections of the waveguides avail- 
Ae ae 7 =o a a7 vi able at the various frequencies did not scale exactly 
FREQ KMC from the 6-kmc case. In addition, for convenience the 
ape poles en OP nso aelG length of ferrite in the various isolators were exactly the 
displacement isolator at 11 kmc. same (5 inches) so as to make use of the same magnet . 


VSWR 
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TABLE. 1 
PERFORMANCE OF FIELD DISPLACEMENT I[SOLATORS 
Maximum Minimum Minimum Ferrite description Applied de 
pe eG forward reverse return magnetic 
( 5 loss loss loss is Magnetization field 
ce (db) (db) (db) Composition (Gauss) (oersteds) 
3700-4200 0.41 1) 2 30 Meg 9CusAl.ssFe1.4Mn.og 1200 590 
5925-6425 0.22 30 30 Meg.sMniFei.s 1960 660 
10, 700-11 , 700 
single slab O25 28 oD Ni.sCu.Zn :Fe;.9Mn.o2 3800 1150 
double slab 0.99 65 Deel 3800 1150 
23 ,500-24 ,470 0.60 34 — Ni.gZn 4Fe;.5Mn_o2 4900 3670 
TABLE Ii an electric field null is obtained at the ferrite face for the 


SCALED AND ACTUAL VALUES OF GEOMETRIC AND MAGNETIC 
PARAMETERS IN ISOLATORS AT VARIOUS FREQUENCY BANDs; 
6-KMC Isorator Is Usep As Basis. QUANTITIES IN 
PARENTHESES ARE SCALED VALUES 


Frequency 5 b h It, 4M, a a 
mc mils mils mils mils | Gauss |*PP J 
oersteds 
5925-6425 180 70 550 1590 1960 660 
3700-4200 RW) WO, V5 [AAO iO, || SOs 
(288) | (112) | (880) | (2544) | 4225) | (412) 
10, 700-11 , 700 IO.) 2a csp 303, MOO, | S800, |) Tso, 
(100) (39) (305) | (883) | (3528) | (1188) 
23 ,500-24 ,470 45, 19, 137, 420, | 4900, | 3670, 
(45) (17) (137) | (397) | (7840) | (2640) 


assembly in all cases. Naturally, as a result of starting 
out with these two nonscaled quantities, we had to 
compensate the other geometric and magnetic quanti- 
ties so as to bring performance at all frequency bands 
into line. This partly accounts for the discrepancies be- 
tween scaled and actual operating conditions shown in 
Table II. 

In addition, the bandwidth over which the 4-kmc 
isolator was tested was considerably greater (on a 
scaled basis) than that for the 6-kmc isolator. Thus to 
obtain good performance over the relatively greater 
band, the geometric quantities had to be altered some- 
what from their scaled values. 

In the 11-kmc isolator, the best ferrite available at the 
time was one with a saturation magnetization approxi- 
mately 300 Gauss away from the scaled value. We see 
from Table II that the position of the slab, 6, was quite 
different from the scaled value to compensate for this. 

The largest 47M, ferrite available at 24 kmc was 
4900 Gauss, considerably away from the 7840 Gauss re- 
quired by scaling. While we used approximate scaled 
values for the transverse geometric parameters, we note 
the applied magnetic field is quite far from the required 
scaled value. Indeed, Table I shows the forward loss for 
this isolator to be quite out of line with forward losses 
obtained in the other cases. We may understand this on 
the basis of the theory of the electric field null developed® 
for the case of a full height ferrite. There it is shown that 


forward direction of propagation if the geometric and 
magnetic quantities satisfy 


T 
oe (uw? — Rk?) tan Rnb 
+ tan k(L — b — 8) = 0 


URm — RB tan ky, 


where all quantities are defined in the Appendix. Clearly 
if the geometric quantities are scaled in going from 
6 kmc to 24 kmc (i.e., Ra, Rm, B—fka, fm, f8 and L, 6, 
6—(L/f), b/f, 6/f, f being the scaling factor), while the 
magnetic quantities are not, so that wand k are different 
at 24 kmc than at 6 kmc, then the equation is not pre- 
served and the electric field null is destroyed. Hence the 
higher forward loss in the 24-kmc isolator. 

Finally, the applied fields would not be expected to 
scale if the internal anisotropy fields are different from 
their scaled values in the various ferrites used. While we 
have no data on the anisotropy fields in the various 
ferrites employed, it is unlikely that they scale according 
to frequency. 

We have made use of scaling only as a guide in ob- 
taining isolators at other frequency bands. It is highly 
likely, however, that if all geometric and magnetic 
parameters were perfectly scaled, comparable perform- 
ance could be obtained in each scaled frequency band. 


APPENDIX 


Scaling for a rectangular waveguide containing a full 
height ferrite slab may be seen from the following con- 
siderations: the phase constant B satisfies the transcen- 
dental equation? 


1/2 B? 
>( + — bn!) cos k,(L — 6 — 26) 
DING pe 6? 
8h, 
Raptlliche vs ri Pe 5 
p@ 
1 (Re? B? 
+>( — —+ kp?) cos ka(L — 8) 
DON pe 6? 


a Re; 


+ cot (Rm) sin k,(L — 6) = 0 


p 
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where 
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Suppose the loaded waveguide produces a phase con- 
stant B at the frequency w. We assert that the phase 
constant at a new frequency fw (f any number) is fB 
provided the geometric and magnetic quantities are 
changed from their old values 6, 6, L, 47M,, Ho to the 
new values b/f, 5/f, L/f, f4mrM,., and fHo at the new 
frequency fw. For then, u, k, p, and @ remain unchanged 
while kg and km become fka and fkm provided 6 becomes 
f® at the new frequency. Indeed, with these changes, 
f@ satisfies the transcendental equation at the new fre- 
quency, as is readily observed. A typical component of 
field is given by? 


Bivia jRa ; 
1 — ——}sin k,(L — 6 — 8) + —cos k,(L — b — 4) |e- tn? 
mp 


BN Ra 
1 — —— }sin k,b — 7 —— cos kab 


m 


4arM, = saturation magnetization of ferrite 
y = gyromagnetic ratio of electron 
w = angular frequency of wave 
woo = yH 
H) = internal dc magnetic field 
k,? = Ww €q {Lo a B? 


w€ML0 


Rane aa ta Ste 
p 


€0, Mo = dielectric constant and permeability of free space 


e = dielectric constant of ferrite. 


Rinp 


and thus the fields are unchanged in value and distribu- 
tion at the new frequency, phase constant, and operat- 
ing conditions provided the geometry x, y is scaled 
down (or up) to x/f, y/f. 
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A Method of Producing Broad-Band Circular Polarization 
Employing an Anisotropic Dielectric’ 
H. S. KIRSCHBAUM} ann S. CHENT 


Summary—A procedure is described whereby it is possible to de- 
sign circular polarizers for both waveguides and in window form to be 
used over a broad band of frequencies. The difference in phase 
constants for two mutually orthogonal E fields while propagating 
in an anisotropic dielectric is combined with the effect due to guide 
wall spacing to obtain a reasonably constant differential phase 
constant for the two fields over a broad frequency band. By properly 
choosing the length of the anisotropic dielectric in the direction of 
propagation, and orienting this dielectric properly with respect to an 
incident linearily-polarized wave, the transmitted wave is circularly 
polarized over a correspondingly broad band of frequencies. 


INTRODUCTION 


ARIOUS methods!“ have been proposed for ob- 

\ taining circular polarization of an electromag- 

netic wave in a waveguide or in free space. Some 
of these methods are frequency sensitive while others 
are not. This paper describes a method for obtaining 
circular polarization over a broad band of frequencies 
by making use of an artificial anisotropic dielectric. 

A wave in a waveguide may be said to be circularly 
polarized if the following conditions are met. The wave 
shall consist of two equal components in space quadra- 
ture (each being a dominant mode) traveling in the 
same direction. The electric vectors of the two com- 
ponents shall be 90° out of phase with each other. The 
equality of the two electric field vectors is necessitated 
by the fact that a circularly-polarized wave can be 
propagated in a waveguide without loss of circularity 
only if the two components that go to make it up are 
propagated with the same velocity. In other words, 
what is required is a guide that is doubly symmetric in 
cross section such as either a square or a circle. 

The method described here makes use of a section 
containing an artificial anisotropic dielectric in the 
waveguide (which is not necessarily doubly symmetric), 
disposed in such a way as to produce circular polariza- 
tion at the interface furthest from the source of energy. 


* Manuscript received by the PGMTT, December 7, 1956. The 
research reported was sponsored in part by the Signal Corps Eng. 
Labs., U. S. Army Signal Corp., through a contract with the Ohio 
State University Research Foundation. ; ’ i 

+ Antenna Lab., Dept. of Elec. Eng., The Ohio State University, 


Columbus, Ohio. ‘ es 
1 “Preliminary Instruction Book for the M4902 Airborne ‘S’ Band 


Circularly-Polarized Radiator,” Radio Res. Lab., Harvard Univer- 
sity Rep. No. 411-1B-71; March 12, 1945. , aby sis 
2 “Quarter-Wave Plate for Broad-Band Circular Polarization, 
M.I.T. Rad. Lab. Ser., Rep. No. 769, McGraw-Hill Book Co., Inc., 
New York, N. Y.; January 28, 1946. ; ; 
3 R. M. Brown and A. J. Simmons, “Dielectric Quarter-Wave and 
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If this interface is followed by a perfectly square or 
round guide, the circularly-polarized wave will propa- 
gate unaltered. 

The use of an anisotropic dielectric is further extended 
to produce a polarizing window which can control the 
polarization emanating from a linearly-polarized an- 
tenna. 


ANALYSIS 


Consider a rectangular waveguide as in Fig. 1, which 
is capable of transmitting both the TE, and TE 
modes. This waveguide is filled with an anisotropic di- 
electric having the property that the dielectric constant 
for the electric vector in the y direction is greater than 
the dielectric constant for the electric vector in the x 
direction. If it is assumed that only the TE;o and TEo: 
modes are transmitted, then the x-polarized wave will 
“see” only e, and the y-polarized wave will “see” only e,. 


x 


Fig. 1—Coordinate system and waveguide dimensions. 


These two waves will behave as though the guide were 
filled with an isotropic dielectric having dielectric con- 
stants e, and €,, respectively. For these two polariza- 
tions, the propagation constants are, respectively, 


B=. Bo / i (=) (1) 
p= po4/ by — (75) @ 


where By =w~/uoeo, kz and k, are the relative dielectric 
constants of the anisotropic medium in the directions 
indicated by the subscript, and ) is the free-space wave- 
length of the wave. 

For the case where €, >e€, and 6 >a, the above relations 
plot as shown in Fig. 2. At some frequency wo the slopes 
of the two curves are equal, and this frequency can be 
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By 


Wey 


Wex Wo 
Fig. 2—Phase constant vs frequency. 
taken as the center of a band of frequencies over which 
(8, —B82) remains very nearly constant. Upon solving for 


the frequency at which the slopes of the two propaga- 
tion functions are equal, there is obtained 


a Les Re 
Ne ee / se (3) 
V boeo( Ry — Rx) os ba ky 


From (3) it follows that 
ht eee 


2 ak, 0b ky 


where Xo is the free-space wavelength corresponding to 
wo. If Xo is fixed at some particular value, then a curve 
such as Fig. 3 can be drawn showing the relationship 
between a and 0 required to satisfy (4). 


Q 


k 
= ¥ 
sae Ve. 


SS 


= 


CEE Xora 
2 Vi, 


Fig. 3—Usable range of guide dimensions. 


It should be noted that this curve is asymptotic to 


Xo V ky/ Re 

2 VWky — ky 
as b becomes infinitely large. The case for 6=infinity 
will be taken up in greater detail later when the subject 
of a polarizing window is discussed. 


Eq. (4) is good provided neither a nor 6 becomes so 
small as to result in the guide’s being cut off at the 
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wavelength Xo. Since a practical design requires b>a, 
and since €,>€z, the guide will cut off first for the 
x-polarized wave. The cutoff wavelength for the «-po- 
larized wave is 


Nos = anh, (5) 


from which it can be seen that the dimension @ cannot 
be less than \o/2/k,. Hence in Fig. 3, only the solid por- 
tion of the a vs b curve can be used. From the point of 
view of the simplicity of the transition pieces from rec- 
tangular to round or square guide, it would be best to 
choose a square cross section for the quarter-wave plate. 
However, with the k, and k, available, such a choice 
would either result in a quarter-wave plate of large cross 
section or it would put wo too close to the cutoff region, 
To permit wo to be somewhere near the center of the 
band of frequencies to be considered, it is usually neces- 
sary to choose a point on the curve of Fig. 3 such that 
at the lowest frequency to be considered the quarter- 
wave plate is not cut off. 

Having established the dimensions of the cross sec- 
tion, a curve of AB vs frequency can be computed for the 
frequency band under consideration. From the average 
value of AG over this band, a length can be chosen for 
the quarter-wave plate such that (AB)ave times the 
length J equals 7/2 radians. 

Up to this point, no consideration has been given to 
the problem of internal reflections within the quarter- 
wave plate due to mismatches at the two interfaces. If 
it is assumed that the guide walls are continued un- 
broken on either side of the quarter-wave plate, as 
would probably be the case for a round or square guide, 
then the complex transmission coefficients through the 
quarter-wave plate for the x and y polarizations of the 
wave are 


: 1 
kiz = 6 
ba + j—( -—) in Bul ) 
cos 8; sin B, 
SADE 
5 1 
Rey = ) (7) 


Le/ Zon Lag 
cos Byl j sin 
By ui - z C : e) Ee 
where #,, and hwy are the complex transmission coef- 
ficients for the x- and y-polarized waves. The quantities 
Zoz and Zo, are the guide impedances for the air-filled 
guide and Z,, and Z, are the guide impedances in the 
region filled with the anisotropic dielectric. From (6) 
and (7) the actual differential phase shift between the 
incident waves and emergent waves will be 


Ad = (by — dz) = arctan ' : E {= | tan B i 
oe PH Aa eae a : 


t | Ee Ba a] B i (8) 
AVRO, Weahiayoras ye 
1 Lies Ziaie 
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The ratio of the absolute values of k,, and wy when 
averaged over the frequency band under consideration 
will serve to determine the relative magnitudes of the 
incident x- and y-polarized waves necessary to produce 
circular polarization in the emerging wave. The incident 
waves must meet the following requirement in order to 
produce circular polarization, namely, 


| I Bi 
Loe 


_ | be 


£3 Re 


(9) 


ave 


where the average is taken over the frequency band un- 
der consideration. 

This requirement follows from the fact that the trans- 
mitted waves must be equal in magnitude in order to 
have circular polarization. 

This same anisotropic material when placed in a 
round guide leads to a problem, the exact solution of 
which is intractable. However, as shown later under the 
experimental work, such a round guide behaves qualita- 
tively in the same manner as a square guide having the 
same cutoff wavelength in the dominant mode. 


THE POLARIZING WINDOW 


As stated above, when 0 is allowed to approach in- 
finity the dimension a@ approaches 


ho Why/ ke 

De ie 
and for the y-polarized wave, the guide would become a 
parallel-plate transmission line. For the x-polarized 
wave, the guide would be infinitely deep. If a number of 
these guides were stacked above each other, there would 
result a configuration as shown in Fig. 4. This could be 
made to serve as a polarizing window when placed in 
front of a linearly-polarized antenna. In Fig. 2, the cut- 
off frequency for the y-polarized wave would become 
zero and the phase constant through the material would 
become a straight line function of frequency. As in the 
case of the quarter-wave plate in the waveguide, (6), 
(7), and (8) apply with Zo, and Zoy equal to the imped- 
ance of free space. The impedances Z;, and Zy are, re- 
spectively, 


V 10/0 
(3) 
2a 
1 Ho 
ene pis 11 
og V ky €0 \ 


Tue ANISOTROPIC DIELECTRIC 


The anisotropic dielectric used in the quarter-wave 
plates described in this paper is an artificial one com- 
posed of alternate sheets of polystyrene and air (or poly- 
foam). The sheets of polystyrene and air are of equal 
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Conducting 
Sheets 


,! 
Fig. 4—Construction details of polarizing window. 


thickness as shown in Fig. 5. If the wavelength used in 
the quarter-wave plates is very much greater than the 
thickness of the sheets of the anisotropic dielectric, the 
dielectric constants of the anisotropic dielectric can be 
computed on the basis of a dc field. Upon doing this it is 
found that 


y 

Ca = ies (12) 
Ep + €f 
€p te 

a : af (13) 


where e, and ¢, are the dielectric constants for x- and 
-polarized fields, respectively, and €, and e; are the di- 
electric constants of polystyrene and air (or polyfoam) 
respectively. These equations are obviously applicable 
to relative dielectric constants also. Tests on the quar- 
ter-wave window placed in front of a horn antenna up 
to frequencies where the thickness of the polystyrene 
was about 4 of a wavelength indicate that the use of 
these dc dielectric constants is justifiable. 


Polyfoam 


Fig. 5—Artificial anisotropic dielectric. 


For a polystyrene and air anisotropic dielectric, the 
e, and ¢€y are 1.438 €) and 1.780 €, respectively. If poly- 
foam is used instead of air, they are 1.470 €) and 1.795 €, 
respectively. 


EXPERIMENTAL RESULTS 


Quarter-Wave Plates 


A number of quarter-wave plates embodying the prin- 
ciples stated above were built and tested. The first of 
these, shown in the photographs of Fig. 6, had the di- 
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(a) 


(b) 
Fig. 6—Photographs of quarter-wave plate. 


mensions listed in Table I, below. This quarter-wave 
plate was designed for a rectangular guide so that ad- 
vantage could be taken of the differential phase shift 
inherent in a rectangular guide regardless of the pres- 
ence of an anisotropic medium. 


TABLE I 
Dimension a ==) SOC 
Dimension } =2.80 cm 
Length / =3.30 cm 
Dielectric thickness =0.159 cm 
Air thickness =0.159 cm 
€&y =1.780 € 
ee = eASoles 


This quarter-wave plate was tested for axial ratio vs 
frequency using the circuit of Fig. 7. The results ob- 
tained are shown in Fig. 8, the anomalous behavior at 
11.1 and 12.45 kmc being caused by resonance within 
the transition pieces. The theoretical performance is 
drawn on the figure as a comparison. 

In order to eliminate the transition pieces, a quarter- 
wave plate, employing the same artificial anisotropic 
dielectric, was built into a round guide of 15/16-inch 
nominal inside diameter. This plate was 4.83 cm long, 
this length being based on the calculated performance 
of a hypothetical square guide having the same cutoff 
wavelength as the actual round guide used. In actuality 
an error in calculation resulted in the length of 4.83 cm 
being about 10 per cent longer than it should have been. 
The axial ratio vs frequency obtained for this quarter- 
wave plate is shown in Fig. 9. In this case, excellent re- 
sults were obtained over a 1.56 to 1 frequency band. The 
theoretical performance of the hypothetical square 
guide quarter-wave plate of length 4.37 cm is shown for 
comparison. These measurements show that less differ- 
ential phase shift per cm of dielectric is obtainable in the 


IRE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 


July 


Oscillator Transition 
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This Can Be Rotated 


This Can Be Rotated 


Rotating Rotating haha Transition Tunable 

Joint Joint Odras Round To Detector 
Rect. Mount 
Fig. 7—Test equipment used to measure axial ratio of 
quarter-wave plates in a guide. 
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Fig. 8—Axial ratio vs frequency for \/4 plate in 2.82.3-cm guide. 
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Fig. 9—Axial ratio vs frequency for \/4 plate in 
15/16-inch circular guide. 


round guide than in the square guide, and that the mis- 
match at the interfaces does not seem to be as poor as in 
the case of the square guide. This may be explained 
qualitatively as follows. Consider the case of both guides 
transmitting waves polarized along the axis of greatest 
dielectric constant (the y axis in these reports). In the 
case of the square guide, there are no field components in 
the x direction and so the wave “sees” only ¢,. In the 
case of the round guide, off the diameter, there are com- 
ponents of electric field in both the x and y directions, 
the field in the y direction predominating. In any event, 
for the purposes of determining the propagation con- 
stant for this wave, the wave “sees” a dielectric constant 
somewhat less than e,. In a similar manner an x-oriented 
wave propagating in the circular guide will “see” a di- 
electric constant somewhat greater than e,. Hence the 
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differential phase shift between x- and y- oriented waves 
will not be as great as it would in the perfectly square 
guide. On the other hand, this effect should improve the 
variation of the ratio of transmission coefficients, and 
result in a better axial ratio over the frequency band 
than would be obtainable in square guide. 


POLARIZING WINDOWS 


In Fig. 10 is shown a photograph of a polarizing 
_ window built in accordance with the above theory. The 
dimensions of the window are given in Table II. Meas- 
urements were made on this window by illuminating it 
with a 6X6 inch horn in which were placed several 700- 
ohm absorbing cards perpendicular to the incident 
polarization. This was done to absorb any cross-polar- 
ized reflections from the window. The incident polariza- 
tion is inclined at an angle of 42.5° with respect to the 
metal strips of the window. This angle was determined 
by test and deviates from 45° because the average trans- 
mission coefficients for the x- and y- polarized waves are 
not equal to each other over the frequency band under 
consideration. The beam from the window was directed 
into a 7-foot cube “dark room” in which was placed a 
polarization-sensitive detector. Measurements were 
made, directly on the beam axis, of axial ratio over a 
band of frequencies extending from 8.2 kmc to 17.2 kmce. 
The results plotted in Fig. 11 show that the axial ratio 
does not exceed 1.5 over a 2.1 to 1 frequency band. It 
further appears that the useful bandwidth of the window 
has not been realized. 


CONCLUSION 


By means of the foregoing design procedure, it is pos- 
sible to build circular polarizers for both waveguides and 
in the form of windows which will operate over a very 
broad band of frequencies. The design procedure did not 
take into account either impedance matching or “end 
effects” in the case of the window. If both of these were 
taken into account either analytically or experimentally, 
it should be possible to secure even better results than 
those obtained here. 


Fig. 10—Photograph of polarizing window. 


TABLE II 


DIMENSIONS OF POLARIZING WINDOW 


AXIAL RATIO 


.64 cm éy=1.79 A. =2.71 cm 
.60 cm €,=1.47 fo=11. kme 
(AB)o 1=82.5° 


THEORETICAL | 


Pea 


20 24 «28 (132 
FREQUENCY (KMC) 


Fig. 11—Axial ratio vs frequency for polarizing window. 
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Errors in a Magic-Tee Phase Changetr* 


RAYMOND M. VAILLANCOURTT 


Summary—This paper recalls the basic properties of a magic-tee 
and how it can be used as a linear phase changer. An analysis of the 
symmetrical magic-tee phase changer is made, which shows that non- 
linearities of the phase shift and amplitude modulation are second 
and higher order effects caused by small mismatches of the structure. 
Also, some qualitative comments are made on the errors of an asym- 
metrical phase changer. Measurements on a phase changer as- 
sembled from ordinary laboratory equipment show that the phase 
shift is linear to better than 1°. 


INTRODUCTION 
Alp HE PRINCIPLE of operation of the ideal magic- 


tee phase changer! is readily understood from a 

consideration of the properties of a magic-tee. 
Many microwave structures have the properties of a 
magic-tee?*® but for the sake of simplicity, the most fa- 
miliar device, z.e., the side outlet “EH” and “H” plane tee 
is used to visualize the principle of operation. The side 
outlet tee shown in Fig. 1, which is assumed to be sym- 


ay 


ee 
by 


Fig. 1—Side outlet “EZ” and “H” plane arm tee. 


metrical with respect to the plane which contains the 
axis of arms 3 and 4, is called a magic-tee when the 
structure is completely matched, 7.e., when a matched 
load is seen looking into any one of the four arms if the 
other are terminated in reflectionless loads. For this 
case, the total power a generator feeds into any one arm 


* Manuscript received by the PGMTT, December 17, 1956. 

+ Canadian Armament Res. and Dev. Establishment, Valcartier, 
Quebec, Canada. 

1G, Saxton and C. W. Miller, “Magic-tee waveguide junction,” 
Wireless Eng., vol. 25, pp. 138-147; May, 1948. 

21m Gre Montgomery, Rewrle Dicke, and E. M. Purcell, “Principles 

of Microwave Circuits, M.I.T. Rad. Lab. Ser., McGraw-Hill Book 
Co., Inc., New York, N. We~ voles) chal2: 1948. 

3 W. K. Kahn, “E-plane forked hybrid-T junction,” IRE Trans., 
vol. MTT-3, pp. 52-58; December, 1955. 


couples in equal amounts into the two adjacent arms 
and not into the opposite arm. If a generator is con- 
nected to arm 4, the coupled waves into arms 1 and 2 
have the same phase in planes equidistant from the sym- 
metry plane, but are 7 radians out of phase if the genera- 
tor feeds arm 3. Now, if two variable short circuits are 
inserted into arms 1 and 2 at a distance Lj, and Ly» 
=I,+(dg/4) respectively from the symmetry plane, all 
the power going into arm 4 will appear at the output of 
arm 3. The wave launched at input 4 splits into two 
equal and in phase waves in arms 1 and 2. Because 
arm 2 is \g/4 longer than arm 1, these two waves are 
reflected back in the symmetry plane with a phase dif- 
erence of 7 radians and will cancel each other in arm 4 
but reinforce each other in arm 3. If the lengths of arms 
1 and 2 are both increased by AL, the total path length 
of the structure has changed by 2AZ, and the phase of 
the output wave in arm 3 has changed by (2AL/\g)2r. 

It should be added that the variable short circuits can 
also be inserted in arms 3 and 4 instead of arms 1 and 2 
and will give a similar relative phase change of the out- 
put in arm 2 when arm 1 is connected to a generator. 

Ten independent equations are available to fix the 
relationship between the 16 parameters required to de- 
scribe the behavior of an unmatched and asymmetri- 
cal hybrid tee. The relations between the various pa- 
rameters are not linear, and solutions for the dependent 
parameters in terms of a suitably chosen set of six inde- 
pendent parameters are practically untractable. On the 
other hand, for an unmatched symmetrical tee, only 
nine parameters are required, and three of them can be 
chosen independently. 

The resulting simplification is such that an analysis 
of the errors in the linearity of the phase change caused 
by reflections is possible. Thus the main part of this pa- 
per is concerned with the analysis of the errors in a sym- 
metrical hybrid tee phase changer. Some qualitative 
remarks on the effects of small asymmetries are made 
at the end of the paper. Experimental results on a typi- 
cal phase changer assembled from commercial com- 
ponents are also included. 


SCATTERING MATRIX OF SIDE-OUTLET TEEF4 


The reference planes in each arm (see Fig. 1) are 
chosen far enough from the junction to be out of the dif- 
fraction field. It is assumed that only the fundamental 
Ho mode can be propagated in the waveguides and that 
the guides have the same uniform cross section and 
negligible losses. 


*C. G. Montgomery, “Technique of Microwave Measurements,” 
McGraw-Hill Book Co., Inc., New York, N. Y., sec. 9.2; 1947. 
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The amplitudes of emergent and incident waves are 
related as follows: 


by ide a1. Sag a ay 
be < Sor Soo S23 So4 a2 (1) 
bs S31 S32 Sa3 Sa4 a3. 
bs 4) Se S30 USay oil 
a, and 6, are so normalized that | 2,.| * avin | by” ate re= 


_ spectively proportional to the incident and emergent 
power at terminal z. Since the structure is assumed 
lossless, its scattering matrix has the properties of a 
unitary matrix, 7.e., 


il 1Or a Ss 


= D, 
0) 1i@® ose & (2) 


4 
ay, SipSsi 
j=1 
The application of these conditions to the scattering 
matrix (1) together with the reciprocity conditions 


Si Se Sa a, oy = 1, yp. 3, 4 (3) 
gives 10 independent equations. But in general, Sn» 
= Umn+jVmn, and hence there are 20 unknown pa- 
rameters. Without loss in generality, the four reference 
planes can always be chosen so that the diagonal ele- 
ments Sp» are pure real numbers.® This choice of refer- 
ence planes reduces the number of unknown parameters 
to 16. This implies that at least six unknowns must be 
measured to evaluate all the others. The above proper- 
ties of the scattering matrix will be used in the analysis 
of the magic-tee phase changer. 


ANALYSIS OF Macic-TEE PHASE CHANGER 


Consider a tee with variable plungers inserted in arms 
1 and 2 at 6; and @ electrical degrees from their respec- 
tive reference planes. If a generator is connected to arm 
4 with arm 3 terminated in a matched load, it follows 
that 


Ca bye 2791 
a= — boe 2782 
Obey = 0 (4) 
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where 
A = (1 + Syre-2#*1) (1 + Syge242) — Syp2eW24 (rt), (6) 


Any actual tee should have very little asymmetry in 
order to have a simple calibration law when used as a 
phase changer. Assuming perfect symmetry makes the 
analysis possible and yet gives a good theoretical insight 
to the actual performance of a good magic-tee phase 


changer. 
From symmetry considerations it follows that 
So3 = — S13 So = Su 
Seog = Si S34 = 0. (7) 


S34 must be equated to zero since, for the M1. mode of 
propagation, the electric field distribution in the cross 
sections of arms 3 and 4 are asymmetrical and sym- 
metrical respectively and hence there is no coupling be- 
tween these arms. Substitution of (7) in (5) and (6) 
gives 

Sat Su U CS rae) 


tb Sya(e-2781 e272) 4. (Sy? — Syq) e724 (0+ 02) 


bs - (8) 
The phase and amplitude of 63/a, can be evaluated from 
(8), but simpler expressions follow without loss in gen- 
erality if the diagonal elements of the scattering matrix 
are made pure real numbers. The Appendix gives the 
derivation of | Sis| ; | Sis| , and | Sia| as functions of Uj, 
U33, and Uy, 1.e., 


told pile aelsan pad ia Weare Ur 
| Sis |2 = ———— ; | Sia |? = —— 
y) 2 
rege hana 
AU 
Wage? i Claas 
Vie? = 9 i. Uy)? aE U2”. 


Substituting these values in (8) the phase of 63/a4 is 
then expressible as 


bz T : 
gad — (ata 7) -y (9) 


a4 


where 


[Uir(sin 20, + sin 262) + a sin 2(0; + 62) + B cos 2(0,; + 62) | 


y = tan’! — 


and the expression for 63; as a function of az is 
leh Gin E a ~ 1 Sue Sie 211 = SaaS ace et? 
+ (—Sy2Si3Se4 — SieSi4Se3 + SisSi4S22 
+ S35 e465 11) e722 1+ 82) | (5) 


5 Montgomery, Dicke, and Purcell, op. cit., p. 149. 


[1 + Uis(cos 26, + cos 262) + a cos 2(6, + 62) — B sin 2(A,; + 42) | 


with 
a= Uy? = Ui + Vise? 
and 
B= 2Uj2Vi2. 


The constant phase angle ¢o is the phase of Si3Si4. The 
modulus of 63/a4 is given by 


bs 2|Sis| | Sia] | sin (61 — 62) | 


a4 ; (1 — a)? + B? + 4U 1, COS (0, = Ao) COS (0, + Ao) 
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(10) 


+ 4| [cos (0, + 62) + Uy, cos (0; — 6) |[a cos (0; + 62) — B sin (0; + 62) |] }3/2 


Substitution of 6.=0,+(7/2)+#, (¢ is a small angular 
error which may arise in the initial setting of the plung- 
ers in arms 1 and 2) in (9) and (10) gives 


b: 
A sme hh 7 (28, ode Vi > (11) 


Since a? and 6? are small quantities of fourth order, the 
square root is of second order and hence the modulation 
amplitude \/a?+? of the output is of second order. The 
output varies at twice the rate of the phase shift of the 
principal wave. For the particular values of Un, Uss, 
Us, used above the modulation of the output is 1 per 


a4 
where cent. 
aeeetes Uy,[sin 26, — sin 2(81 + t) | + asin 2(26, + £) + B cos 2(26 + 2) 
Be frees 0, Gos 2 (6) to) aw Gos 220 2 LB see 
bs 
a 


2| Sis| | Sia | cos ¢ 


(12 


The third term, yi, of (11) is a contribution of second 
order provided t, Ui, Ui, and Viz are small. Moreover, 
if ¢ is a small error, its variation 6¢ with plunger setting 
will be of second order. Thus to a second order of ap- 
proximation, the phase of 63/a4 is a linear function of 6; 
since @o is constant. The extrema of the third term 
should be very near and of the order of magnitude of 
those obtained for the case t=0. For f=0 


a sin 40; + B cos 46; 


ee 1 — acos 46; + Bsin 46, ve 
and y; has extrema for 
1 1 aie 
0, = er cos—! fo ae alse os | \ ‘ (14) 


To give an order of magnitude of the expected error due 
to di, let Uy = U33= Uy=0.1 (this implies a residual 
vswr = 1.22) and ¢; max~0.57°. When ¢=0 

2 S13 | | S14 | 
{1 + a? + B? — 2a cos 46, + 26 sin 40,} 1/2 


bs 


a4 


1 Pate DAR 
2| Sis| | Sis {1-5 t +0) + VEEP 


=> co a Seas sin 10, | -b sees \ 
Va? + B? A/a GF 
= 2| Sis| Sis| [1 + Va? + 6? cos (40: + 1) (15) 


where 


S 46, 


¢; = tan! B/a. 


{(4 — a)? + 62+ 4Uy, sin ¢sin (20, + t) + 4[sin (26, + 4) + Uysin t][a sin (26, + #) + B cos (26; + f) jue 


EXPERIMENTAL VERIFICATION 


To use a side outlet tee as a phase changer it is neces- 
sary to match it reasonably well (vswr<1.3) and to 
have the plungers in the “through” arms respectively at 
6, and 6:+90 electrical degrees from the symmetry 
plane. 

For narrow-band operation, traveling screw tuners 
may be used in the “EZ” and “H” plane arms. A residual 
vswr~™1.03 is readily obtained provided the load termi- 
nations in the other three arms have small “phasable” 
reflection coefficients. Long tapered wood loads which 
can be slid in the waveguides do very well. 

To fix the respective locations of the plungers in the 
“through” arms, one reflectionless load termination of 
the through arms is replaced by a plunger and the posi- 
tion of one minimum of the electric field in the input 
slotted line is measured. Repeating the same experiment 
for the other arm with another plunger gives the differ- 
ence in the electrical length of the two through arms. 
From this result and the measure of the guide wave- 
length (not the slotted guide wavelength), one of the 
two plunger arms can be made 90 electrical degrees 
longer than the other. 

Because the structure was not perfectly matched, the 
locations of the electric field minima in the slotted line 
were averaged for open and short-circuit conditions at 
the output by moving the output plunger a quarter of a 
guide wavelength away from its initial position for every 
setting of the plungers in the through arms. 

The results of the measurements are shown in Table I. 
The maximum variation of the incremental phase shift 
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Vaillancourt: Errors in a Magic-Tee Phase Changer 


TABLE I 

Din 

cm Vs Ag 
0.000 0.00 
0.254 19.6 19.6 
0.508 39.5 19.9 
0.752 58.8 19.3 
1.016 77.9 19.1 
1.270 97.7 19.8 
1.524 117.5 19.8 
15.753 137.3 19.8 
2.032 156.8 19.5 
2.286 176.4 19.6 
2.540 195.8 19.4 

f=9193 +5 me. 


gAgs =2.327+0.002-cm wavelength in the guide. 
2Agp = 2.340 +0.005-cm wavelength in the slotted line. 
vswr for “E” and “H” arms~ 1.03 
2L;=incremental path length of the through arms. 
¢=relative phase shift in the slotted line. 
A¢=incremental phase shift in the slotted line. 


is within 0.8°. The estimated maximum error in the 
plunger setting Z; is +0.001 cm and in the location of 
the field minimum +0.005 cm. Thus the estimated 
maximum error in ¢ is +0.4°. In view of this limited ac- 
curacy in the measurements it would seem that the 
phase shift was linear to better than 0.8°. 


CONCLUSION 


It has been shown that the symmetrical magic-tee 
makes a lossless linear phase changer whose phase 
linearity is affected to a second order for reflection co- 
efficients equal to or less than 0.1. 

The fact that isolation between the “E” and “H” 
plane arms is not infinite although usually greater than 
40 db (| Ssa| <0.01) introduces additional errors. This 
coupling is caused by a certain amount of asymmetry 
for which S34 may be considered as a measure. Reference 
to (5)—(7) shows that when the tee is not perfectly sym- 
metrical, $3440 and the term (S13S14.S02 + So3.Se45S11 
— Sy.S13S04 —Si2514523) does not completely cancel. This 
term may however be expected to be small for slight 
asymmetry since it does cancel for a symmetrical tee. 
The main effect of small asymmetry is then to introduce 
two additional terms the order of magnitude of which is 
0.02 and varying in phase with respect to the dominant 
term S}3.S,4(e—24 — e—?/) in the numerator of (8). The 
dominant term has a magnitude of approximately unity, 
so the beating of these terms could be expected to pro- 
duce an additional phase error of around + tan! 0.02 
or +1.14° if both perturbing terms add in phase. Thus 
a maximum phase error or departure from linearity of 
2° may occur if all the errors add in phase and are 90° 
out of phase with the dominant wave. The experimental 
results show that in practice errors of less than 1° can be 
obtained using commercially available magic-tees. 
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APPENDIX 

Since [.S] is unitary, it follows that 
OM Sie tree leSag| 4 1 Sua |? (16) 
2| Sis]? + Us? = 1 (17) 
2| Su]? + Us? = 1 (18) 
Uu(Si* + Si) — | Sis? + | Sul? = 0 (19) 
Si3*(Uur — S12) + Si3U33 = 0 (20) 
Si(Ou + Sta) + SuUas = 0. (21) 


These six independent equations have nine unknown 
variables. It is then in general necessary to know the 
values of at least three of them to solve for all the 
others. Because the diagonal elements are the reflection 
coefficients measured in the four reference planes, they 
seem an appropriate choice of values to use. 


From (17) 


1 — Us3;3? 
| Sis |? = ———- (22) 
2 
From (18) 
1 cr U 44" 
| Sia]? = (23) 
2 
From (19) 
U D vas GE 74 
Sie at Si2* = sTaOMR i 
Now 
Sy = Uy + jVie, 
Sie + Sir* = 2U 2, 
and hence 
U Ohaees (iy 2, 
Das B 44 33 d (24) 
4U 1 
Substitution of (22) and (23) in (16) gives 
U 2 U 2. 
| Siz |? =— , a U 41’. (25) 
But 
| Sie [2 — U x2" a Vie? (28) 
whence 
U 2 U 2 
Vi = Cea Ui eon ae (27) 
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Excess Noise in Microwave Crystal Diodes Used as 
Rectifiers and Harmonic Generators” 
J. M. RICHARDSON} anp J. J. FARISt 


Summary—Excess noise produced by microwave excitation of 
silicon crystal diodes was studied for operation of the crystal as a 
detector and as a microwave harmonic generator. The noise appears 
at the detector terminals and also as noise sidebands of the micro- 
wave harmonic, thus degrading the spectral purity of the harmonic 
relative to that of the fundamental. Possible models of the processes 
involved are presented. Difficulties and technique of measurement 
are discussed. Observations for 1N26 crystals, used as detectors, 
doublers, and triplers, and excited by X-band power in the range 8 to 
100 mw are presented, showing limitations on spectral purity set by 
the process of noise production during harmonic generation. 


INTRODUCTION: MOTIVATION AND PROBLEM 


example, microwave spectroscopy and interferome- 

try, it is important to know the spectral purity of 
the sources used because the problem is usually not that 
of the minimum detectable signal but that of the mini- 
mum detectable change in a strong signal. If the trans- 
mission and detection portions of an experiment have 
been made sufficiently noise free, the source noise tem- 
perature may set the limit of precision in the measure- 
ments. 

The most important present source of monochromatic 
waves in the millimeter region is the crystal diode har- 
monic generator. This device is, for convenience, also 
frequently used in producing low-level centimetric 
waves from uhf inputs. Since there were indications 
that the output harmonics were rather noisy, a study of 
the multiplication process and the associated noise was 
undertaken. The existence of excess noise in a crystal 
detector having the familiar 1/f spectral density is al- 
ready known, and it was asked whether these phe- 
nomena were related. For example, the excess noise cur- 
rents flowing as a result of the fundamental microwave 
excitation and its rectification may become modulated 
onto the harmonic by virtue of the nonlinear nature of 
the crystal and appear as noise sidebands with spectral 
density and magnitude characteristic of the low-fre- 
quency excess noise. 

The problem may then be stated as: Given a crystal 
diode n-harmonic generator excited by power P;, at 
angular frequency w, investigate the noise sideband 
power associated with the mth harmonic, S,, if any. S, 
may be investigated as a function of P; and frequency; 
and its relationship to So, the excess noise power of the 


N precision work in microwave measurements, for 


* Manuscript received by the PGMTT, December 19, 1956. Part 
of this work was carried out under a program of research and devel- 
opment in basic instrumentation cooperatively supported by ONR, 
AFOSR, AEC, and NBS. 

+ Natl. Bur. of Standards, Boulder, Colo. 


crystal used as a detector, may also be investigated. The 
nature of .S, may be studied in any convenient way, but 
perhaps the simplest is by demodulation of the spectrum 
in the neighborhood of nw. 


MODELS 


The problem can be considered as a generalized 
mixer problem in which the exciting microwave may be 
thought of as the local oscillator.! Fig. 1 illustrates the 
model and Fig. 2 illustrates the spectral relations in- 
volved. The harmonic generator is fed at the terminals 
w with high level microwave power P; in the range 10 to 
100 milliwatts. Rectified power or harmonics will ap- 
pear at the terminals dc, 2a, - nw. For any particular 
harmonic generator, say a tripler, the terminals at the 
desired harmonic are connected to a matched load and 
an attempt to reactively terminate the other microwave 
terminals is made. The de termination is a resistor pro- 
viding self bias for the crystal. 

In the absence of input power other than at the fre- 
quency w, no other outputs except those just mentioned 
will occur. If, on the other hand, power is fed to the 
terminals 6, perhaps by the internal production of a 
component of excess noise in the crystal, converted 
power will appear at all signal and image terminals, 
nw+.In Fig. 2, the power at 8 and nw +8 is shown dis- 
tributed over a small band to emphasize its possible 
noisy character. The band represents the limiting band 
of the apparatus. Interest is focused on near-carrier 
noise components (8 small) for two reasons: 1) The se- 
verity of the effect is greater for small 6 due to the 1/f 
excess noise spectrum, and 2) it is desired to analyze 
experiments which may involve slow signal modulation 
of the carrier nw, as for example, by 30 cps Stark or Zee- 
man modulation of a microwave spectral line. Side- 
bands due to this signal modulation would then be 
superimposed on the noise-sideband power, to give some 
signal-sideband/noise-sideband power ratio. 

Although the above description, in terms of general- 
ized mixer theory, is undoubtedly a correct and poten- 
tially fruitful model, analysis of the problem was not 
attempted in these terms. This decision was taken be- 
cause of the formidable problem of specifying the ad- 
mittance matrix of the mixer at a particular P, and also 
as a function of Pi, and of specifying the terminating 
impedances of all the terminals. 


1H. C. Torrey and C. A. Whitmer, “Crystal Rectifiers,” Rad. 
nai Ser., vol. 15, McGraw-Hill Book Can Inc., New York, N. Y.; 
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Fig. 1—Generalized mixer showing high- and low-level terminals 
at various frequencies of interest. 


Fig. 2—Spectrum of interest in the detector-harmonic generator prob- 
lem. The lines at 0, w, 2w, - - - , mw are spectral lines whose heights 
have been depicted proportional to the power they represent. 
The small bands represent spectral densities of noise power whose 
integrated values over a specified band are denoted by So, Sz, etc. 


An alternative model was, therefore, used for guid- 
ance. Some current-voltage characteristic, i(e), valid at 
all frequencies of interest,? is assumed, and is expanded 

_about the origin in the power series 


N=) iane™ (1) 
In this series the applied voltage, e, has the form 
e = Ey + FE, cos wi + Bit) (2) 


where Ep is the self bias developed across the diode 
terminals, F; cos wt is the applied microwave voltage at 
angular frequency w of amplitude £; corresponding to 
P,, and B(t) isa time dependent virtual applied voltage 
descriptive of the excess noise generated by the crystal 
upon excitation. B(t) may be thought of as narrow-band 
noise near the frequency 6, whose magnitude is de- 
scribed by B. Ey and 8 are functions of E; (or P,). 
Substituting (2) into (1) and retaining several terms 
in the power series because of the highly nonlinear char- 


2 This assumption is of course not strictly valid, as seen from the 
usual equivalent circuit of the crystal diode. 
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acteristic and large excursion of e, it is possible to write 
expressions for P, and S, in terms of P;, Mo, and B, and 
the coefficients dm. The dependence of Hy on P; is ob- 
servable from the experiment. The result is that 


Sn SGA 
and 
Sale = te (3) 


where G, and H, are functions of Pi, Ho, and the dm. 
Eq. (3) have such form as to permit various dependences 
of the various S, with P. 


PROBLEM OF MEASUREMENT 


Although the source is made as noise free as possible, 
it is unavoidably noisy to some degree, and care must be 
taken that the observed noise is not demodulated source 
noise, S;, (that is, noise associated with P,) in the case 
of the detector experiments, and does not arise from 
source noise modulated onto the harmonics in the case 
of the multiplier experiments. Source noise would appear 
at the detected output of the crystal with relative ampli- 
tude determined by S,/P; and the shape of the P;-Po 
curve. Analysis using the above described models sug- 
gests that noise arising from S; would appear on all har- 
monic outputs in relative strength determined by S;/P1, 
and the shape of the P;-P, curve. 

Thus it was found necessary to carry out the measure- 
ments using techniques that balance out the effects of 
the noise modulation of the source. The rejection ratio 
required to insure that the measurements obtained were 
not influenced by the source noise modulation could be 
estimated by making measurements with the system 
both balanced and unbalanced. 

As a further precaution, demodulation of the har- 
monics was always carried out using a bolometer, an 
excess-noise free device, so that it was certain that the 
bolometer output, if well above its own thermal noise 
level, consisted essentially of noise demodulated from 
the harmonic and its noise sidebands. 


‘TECHNIQUE OF MEASUREMENT 


Although interest was focused on the generation of 
millimeter waves from centimetric waves, measure- 
ments were conducted at X band to enjoy the conven- 
ience of complete instrumentation at the fundamental 
and the first few harmonics. Type 1N26 crystals, a use- 
ful type in millimeter wave generation, were used. 

The source of X-band power used was a klystron, 
type SM X-32, operating as the final stage of a frequency 
multiplying chain. The first stage of this chain was a 
temperature-controlled quartz crystal oscillator. 

Measurements of So, the excess noise in detectors, 
were made with the apparatus of Fig. 3. The detector 
mount terminating the branch produced by the shunt 
tee was used for two purposes. First, a bolometer could 
be placed in this mount to monitor the noise modula- 
tion on the output of the source. Second, a crystal 
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Fig. 3—Microwave circuit used for investigation of 
excess noise in detector diodes. 


placed in this mount was used to obtain a signal to be 
introduced into a differential amplifier with the signal 
from the crystal under test. Thus the input to the dif- 
ferential amplifier consisted of in-phase detected source 
noise from each crystal and incoherent excess noise from 
each crystal. The former noise did not appear in the 
amplifier output because of the differential action, while 
the latter combined according to an rms rule. 

The ferrite modulator was used for two purposes. 
First, it was used to produce modulation at the center 
frequency of the narrow-band amplifier to tune the test 
crystal to eliminate reflections. Second, it was used to 
produce modulation to enable adjusting the attenuator 
preceding the balancing detector for maximum dis- 
crimination against source modulation. Rejection ratios 
of over 100 times were easily obtained for the detected 
voltage due to source modulation. For all data taken, 
this rejection ratio was adequate. 

Because of the balancing technique employed, in 
order to obtain a single measurement of So, the excess 
noise in the crystal under test, three measurements were 
needed. First, the noise present in the balanced system 
described above was measured, then the noise present in 
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Fig. 4—Microwave circuit used for investigation of noise modulation 
introduced by crystal diode harmonic generators. 


the test crystal and in the balancing crystal were each 
measured separately. From these three measurements, 
enough information was obtained to find the excess noise 
generated by the test crystal alone. The crystal under 
test was always operated with a load resistance of 
10,000 ohms. 

Measurements of Sy and $3 and P: and P3 for doublers 
and triplers were made with the apparatus of Fig. 4. 
Again, it was found necessary to eliminate the effects of 
the source modulation by using a differential method. 
The bolometer terminating the branch formed by the 
shunt tee and the bolometer demodulating the multi- 
plied power both produced signals due to the noise on 
the source output, whereas the bolometer detecting the 
multiplier output also detected the noise introduced by 
the process of multiplication. Thus the differential ac- 
tion of the circuit rejected the in-phase noise produced 
by modulation of the source and recorded only the noise 
produced by the multiplication. As in the measurements 
of So, it was necessary to adjust the power demodulated 
at the balancing bolometer to obtain maximum rejection 
of the in-phase signal. The ferrite modulator was used 
to produce a modulation for the purpose of making this 
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adjustment. It was found that a rejection of over 50 
times could easily be obtained. For all of the data taken, 
this ratio was determined to be adequate. 

Finally the multiplied power could be measured with 
a power bridge at the bolometer following the multiplier. 

The multiplying crystals were always operated with 
the bias resistance that produced the maximum multi- 
plied power output. When doubling, this resistance was 
usually in the range from 50 to 200 ohms. For tripling, 
a broad plateau of maximum power occurred for bias 
resistance of about 5000 ohms and above, therefore in 
this case the bias used was always 10,000 ohms. 

The measuring apparatus is shown in Fig. 5. The 
audio-frequency harmonic analyzer was used as a nar- 
row-band af voltmeter, centered at either 270 or 540 
cps. Its output was averaged in a circuit with time con- 
stant of approximately 25 seconds and recorded with a 
recording milliammeter. Each measurement was taken 
by operating the recorder at least three minutes and the 
resulting trace was averaged visually. 
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PREAMPLIFIER} | | pasa 


AVERAGING 
CIRCUIT 


D.C. AMPLIFIER 
RECORDER 


MONITOR 
OSCILLOSCOPE 


> 
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Fig. 5—Block diagram of noise recording apparatus. 


RESULTS 


Results are presented as the square of a modulation 
coefficient because this is a measure of the spectral 
purity of the wave involved, namely the ratio of noise 
sideband power in a given band, S, to the carrier power, 
P. As described above, S is the remaining noise after 
elimination of the effect of source noise. We have 


Mn? = 2Snf Pr (4) 


where n=2, 3 for the doubler and tripler, and ~=0 for 
the detector. The definition of modulation coefficient 
used agrees with Goldman* for a complex waveform, 
and is on the basis of unit bandwidth. Values of S, given 
by (4) may be used to obtain noise temperatures of the 
.multipliers as sources if desired. 

Table I shows the variations of spectral purity with 
crystal, input power, type of operation, and frequency. 
In general, a crystal which is noisy as a detector (large 
my) is also noisy as a multiplier. This result is expected 
if the addition of noise in the multiplication process 
arises from excess noise. The value of m,,” usually in- 
creases with P;, for 7=0, 2, 3 although there are excep- 
tions. This fact substantiates at these high power levels 


3S. Goldman, “Frequency Analysis, Modulation and Noise,” 
McGraw-Hill Book Co., Inc., New York, N. Y.; 1948. 
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TABLE I 

Excess NotsE MopuLATION COEFFICIENTS TOGETHER WITH OuT- 
PUTS FOR 1N26 CrystaLs USED AS DETECTORS, DOUBLERS, 

AND TRIPLERS 


Crys-| P: db 


Lapabielent 10127292 | 10227729? | 101272? | 10127732 Vo P» P3 

No 100mw/i2t270~\at270~ at540~/at270~| (volts) | (uw) | (uw) 

2 1 See, 0.6 0.3 2 4.8 850 | 510 
3 ae 0.6 0.4 2 LW 550 | 320 
5 4.8 0.6 0.4 1.0 oO 330 | 220 
a 4.8 0.6 yal) 

60 1 25.0 0.8 0.6 18.0 ees 830 | 280 
3 13.0 O25 0.5 8.4 2.4 550 | 240 
5 14.0 0.8 0.4 5.8 Dn 340 | 170 
i 13.0 3.6 Bop 100 
9 14.0 Pl 

61 1 200.0 3.6 (ae el AOKO Mex, 980 | 250 
3 94.0 5.8 4.0 23.0 2.4 610 | 200 
S 590 10R2 8.4 27.0 2.4 330 | 120 
7 A0FON 1320 13.0 30.0 2.2 170 59 
9 30.0 | 19.0 10.9 PEAY 74 
11 15.0 | 16.0 1.6 26 

62 1 72.0 20 hoe 32.0 2 <4: 770 | 280 
3 42.0 BEY) 1.4 20.0 Dee: 530 | 220 
5 49.0 1.4 1.0 15.0 2.4 330 | 150 
7 50.0 1.4 0.6 8.4 Ded 170 83 
9 49.0 21 
11 35.0 1.9 

63 1 29.0 0.4 0.3 19.0 yah 700 | 270 
3 26.0 0.3 On2 19.0 2.6 460 | 200 
5 24.0 OZ gil 19.0 2.5 290 | 110 
7 29.0 18.0 2.4 46 
9 30.0 Ded 
11 36.0 2.0 


the strong dependence of excess noise production on 
input power quoted by some authors’ as P,? at lower 
power levels. It also is seen that the tripler output is less 
pure than the doubler output for all crystals and input 
powers. This fact is understandable in terms of the re- 
duced efficiency of their harmonic production, whereas 
excess noise production for a given input power is un- 
altered. 

Table I also shows that the noise modulation at 540 
cps on the doubler output is smaller than that at 270 
cps. This result can also be explained by the assumption 
that excess noise is the primary cause of the multiplier 
noise. It is known that excess noise has a power spec- 
trum that varies approximately as 1/f, so that modula- 
tion sidebands arising from this noise should have the 
same power spectrum by (3). 

Finally, m3? seems to bear a closer relation to mo? than 
to m,” as to both magnitude and dependence on P. 
This behavior is possibly related to the fact that mo and 
mg; were observed using the same bias resistance, where- 
as my was observed with a lower bias resistance. 

In the practical case of a fundamental source of modu- 
lation coefficient m, the total relative noise in the out- 
put would probably be given by adding a term of the 

4M. W. P. Strandberg, H. R. Johnson, and J. R. Eshbach, “Ap- 


paratus for microwave spectroscopy,” Rev. Sci. Instr., vol. 25, p. 776; 
August, 1954. 
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order of m,? (due to an excess noise-free multiplier) to 
m,,” (due to excess noise only). Thus the amount of deg- 
radation of signal purity depends on the characteristics 
of the particular fundamental source used. In this work, 
the source had m,=0.65 X107® (m,?=0.42 X107!") so 
that the degradation was high with all tripler crystals 
and “noisy” doubler crystals, but not severe with 
“quiet” doubler crystals. 

The rather important question of whether the multi- 
plier noise produces fm sidebands in addition to AM 
sidebands was not answered by this study, since only 
the AM sidebands were recovered. However, if as sug- 
gested, the mechanism is that of excess noise, the pro- 
duction of fm sidebands should not occur. 


CONCLUSION 


The results of an investigation into the degradation of 
spectral purity by the processes of detection and har- 
monic generation in crystal diodes have been presented. 
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No completely satisfactory model of the process of har- 
monic generation has been developed, but two possible 
approaches were presented. The measurements are diffi- 
cult because of the need to avoid the influence of residual 
source modulation on the results. An X-band source of 
equivalent noise modulation index of about 0.7 X10 
for unit bandwidth centered at 270 cps from the carrier 
was satisfactorily used. For 1N26 crystals as detectors, 
doublers, and triplers excited with X-band input power 
from 8 to 100 mw, noise dependent on the individual 
crystal, input power, and frequency is added to the out- 
put. The noise modulation introduced by the multipli- 
cation process has characteristics similar to the excess 
noise observed in crystal detectors. 

Further investigation, involving the examination of 
many more crystals, more complete spectral density 
data of the various noise powers, and fm noise would be 
desirable as an aid to the designer of precision experi- 
ments in the microwave and millimeter wave region. 
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Exponential Transmission Lines as Resonators 


and Transformers* 
RABINDRA N. GHOSE 


Summary—An attempt has been made to analyze the theory of 
an exponential transmission line from its complex reflection coeffi- 
-cient’s standpoint and to indicate how the characteristics of an ex- 
ponential line can be completely represented for any frequency with 
the help of the Smith Chart. It is shown that the optimum design 
parameters of an exponential transmission line which may be used as 
a transformer, with a frequency-sensitive load at one end, can be 
determined with the help of the Smith Chart and some derived equa- 
tions. This paper also includes a study of the coaxial type exponential 
line which can be used as a series or parallel resonator. Theoretical 
expressions for the attenuation constant, stored energy, and Q for 
such types of resonator have been derived. Also indicated in this 
paper is the possibility of replacing the uniform-line coaxial-type 
resonators in many microwave and uhf wave filters by the ex- 
_ ponential-line resonators, particularly when a large power-handling 
cpacity is warranted. 


INTRODUCTION 


N RECENT YEARS the exponential transmission 
| line has found wide application in microwave net- 

works as a matching device suitable for matching 
two unequal impedances over a wide band of frequen- 
cies. A general analysis of the exponential transmission 
line has been made by Burrows! and Schelkunoff.? The 
purpose of this paper is to analyze the theory of the 
exponential line from its reflection coefficient’s stand- 
point and to indicate how the characteristics of an ex- 
ponential line can be represented with the help of a 
Smith Chart which is primarily designed for uniform 
transmission lines. The possibility of using the expo- 
nential transmission line in the form of a resonator is 
also discussed in this paper. 


GENERAL EQUATIONS AND THEIR SOLUTION 


For any transmission line system the differential 
equations for the voltage and current can be repre- 
sented by 


aV 
dz 
al 
dz 


+ Z(z)I = | 
(1) 
+ V(2V = 0 


* Manuscript received by the PGMTT, January 30, 1957. Pre- 
sented at the IRE-WESCON Conf., Los Angeles, Calif.; August 21- 
24, 1956. 

+ The Ramo-Wooldridge Corp., Los Angeles, Calif. Formerly 
with Radio Corp. of America, Camden, N. J. 

1C. R. Burrows, “The exponential transmission line,” Bell Sys. 
Tech. J., vol. 17, pp. 555-573; October, 1938. 

2S. A. Schelkunoff, “Electromagnetic Waves,” D. Van Nostrand 
Co., Inc., New York, N. Y., p. 222. 


where 


V =voltage across the transmission line at gz, 

T=current in the transmission line at z, 

Z =equivalent series impedance per unit length of the 
line, 

Y=equivalent shunt admittance per unit length of 
the line. 


For a nonuniform transmission line, Z(z) and Y(g) are 
functions of zg; 7.e., they are functions of the position 
along the line. Eq. (1) will usually give rise to a set of 
second-order differential equations 


@V od dV 
— — [In Z|] — —-ZYV = 
dz? —s dg dz 
(2) 
Clad 
— — [In 'V] — ZYI =0 
scale, ae dz 


Walker and Wax’ have shown that these second-order 
differential equations can be converted into a single 
first-order nonlinear differential equation 


fe ype Ol aaa aG (3) 
dz 2 dz 
where 
jee LZ, ae 
Viiee7, 
and 
y =i /ZY. 
For an exponential transmission line,‘ 
Zig = hee (4) 


where & is the characteristic impedance of the line at 
z=0, and a is the rate at which the characteristic im- 
pedance changes exponentially with the position z. 

From (3) and (4) one obtains the differential equation 
for an exponential line 


dp 


a = 
Party lah oN eile a her y (5) 


3L. R. Walker and N. Wax, “Non-uniform transmission lines and 
reflection coefficients,” J. Appl. Phys., vol. 17, pp. 1043-1045; De- 
cember, 1946. 

4 For an exponential line whose Zo decreases from 2=0, a in (4) 
can be replaced by (—a). 
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If one is interested in the coaxial form of exponential 
line through which TEM waves propagate without cut- 
off, y can be regarded as a constant function of the po- 
sition z, since both a and 6 =27/) are independent of sz. 
This is also true for the parallel-wire transmission line 
and the strip transmission line carrying TEM waves. 

Thus, making use of the following transformation 


S? = (2 f+ 1) (6) 


a 


z - S / 
2 


aSZ : 
tanh —— + tanh SC 
2y Bs 
—+5 = = (8) 
a aSz 
(1 + tanh a tanh sc) 


dy 
Re ple 
a 


one obtains from (5) 


and 


p(z) = — 


where C is the constant of integration which has to be 
evaluated from the boundary condition. 

Let the boundary condition be so assumed that 
p(l) =po, where po is known. From (8), then, 


2y aS] 
po + — psa 


fail SOs i (9) 
2y aS] 
S + @ =e -*) tanh cane 


a 


2 Le 
p(z) aoe aE (= ay Zu), 


where Zn is the input impedance of a uniform transmis- 
sion line with normalized surge impedance of one and 
which is terminated at the load end, z=/, by an im- 


pedance 
= Po av 
7 any 
; G oS 


As the equivalent Zin can readily be obtained from 
the Smith Chart when a, y, and po are known, the effect 
of the exponential-line transformer in terms of the re- 
flection coefficient and vswr at the input end can be de- 
termined without too much laborious computation even 
when a frequency-sensitive load with arbitrarily varying 
reflection coefficient is connected at the load end. It may 
be remarked that while computing Zin, the electrical 
length of the line should be considered as 8/7, where 


and 


(10) 


(11) 


p= g4/1 ae 
2B? 


for a lossless exponential line transformer. 


(12) 


IRE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 


July 


To enable the transformer to work well beyond the 
cutoff frequency, it will be desirable to choose @ such 
that 

Aor 
CERES (13) 
oN 
\ being the wavelength corresponding to the lowest fre- 
quency in the passband which has to be transmitted 
through the transformer. The fictitious impedance Zz 
in (11) is, in general, complex. Assuming the trans- 
former is lossless and is operating well beyond cutoff 


y = 18 
and 
= B po 
Z,=——i=: 14 
l B L 26 ( ) 


In order to make use of the method described, it is es 
sential that Z, should be in the right half-plane, as 
otherwise one cannot make entries on the Smith Chart 
to compute Zin. To ensure this, one can set 


_ Sie pe < exo i(@ = =) It >0 (15) 


where @ is the phase angle of the reflection coefficient po 
at the load end of the transformer. That is, 


a> || acos(o+*) 


But | po| <1 and the value of cos (@+7/2) <1. Hence, 


Tv 
| po | cos (0+ =) <1 


and 
B>— 
2 


But this condition has already been assumed in order to 
enable the transformer to operate well beyond the cut- 
off frequency. Hence no difficulty will be experienced in 
computing p(0) according to the method described 
above. An equivalent circuit describing the method of 
computing p(0) is shown in Fig. 1. It should be recog- 
nized that this method will be of considerable help in 
synthesizing the design parameters of a transformer for 
optimum p(0) over a frequency band with any specific 
load connection.® 


EXPONENTIAL LINE AS A RESONATOR 


From the preceding analysis, it appears that a stand- 
ing wave can be maintained in a section of lossless ex- 
ponential line when po is chosen +1. This suggests the 
possibility of the use of the exponential line as a reso- 
nator. 

5 Schelkunoff, op. cit. The reflection chart shown in Fig. 7.11 
can be used to determine the complex reflection coefficient when the 
normalized impedance of the arbitrary load is known. Similarly, the 


input impedance of the transformer for an arbitrary load can be 
determined readily from Fig. 1 above and the reflection chart. 


S =) 
as is 
Rg 
wn 
NI 
is} 
mer) 
< fll ees 
a COU. 
a) 
2 fede 
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foo) 
+ 
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als 
p(0) = -i5F+Sz, 


Fig. 1—Equivalent circuit of an exponential transformer. 


Let it be assumed that one is interested in using a sec- 
tion of an exponential line as an infinite-impedance or 
parallel-resonance type resonator when one end of the 
line is short-circuited at z=/, as shown in Fig. 2. This 
is analogous to a \/4 cavity for a uniform line 


(16) 


20 eae iP 


Zea) 


Fig. 2—Exponential resonator with Zp increasing with z>0. 


From (10), 


aS 


| a 


where k is any integer. The required length for parallel 
resonance is 
(2k — 1) 


CES 


and 


r 


a? 1/2 
(1 = <, 
4p? 


where a lossless line is assumed, such that y =7@. Simi- 
larly, if one is interested in using the exponential line as 
a series resonator, 


(17) 


Vronin oF 


nr (18) 


aS 
Po | z=1 aie 7 p(0) ee 1, 9 ae 


where 7 is any integer. 
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The length of the resonator is 
ATH 
bs = ae eae 
(4B? —_ a) 112 
and 
oN 
Lente =< (19) 


me \ 1? 
2 (1 3 =) 
48? 


Fig. 3 shows the required length of the resonator for 
series- and parallel-type resonance. 


0.50 1.0 
0.40 0.8 
A ee 
Fe 0030 0.6 § 
E E 
¥ an 
am oF 
0.20 /_—__— —|— 0.4 
0.10 a 0.2 
1.0 


x/B —> 


Fig. 3—Minimum resonant length of series and 
parallel-type exponential resonators. 


STORED ENERGY, WALL LOSSES, AND Q OF 
AN EXPONENTIAL RESONATOR 


From the differential equation of voltage in the ex- 
ponential transmission line whose characteristic imped- 
ance increases with z>0, one obtains 


AZ 2s AZ 
V(z) = A exp (= — iB) ae exp (= =" ite) (20) 


where A and B are arbitrary constants. When a shorting 
plate is placed at z=], 


and 


(21) 


216 


Let the current at any z be represented as 
I(p, 2) = If(z). 
For the TEM mode, 


nl f(z) 
2p 


b(z) Tf Z) 
V = ih mie) dp 
a(z) 2mp 


sin (Bl — Bz) 
(2) J” 
[a 
a(z) 
n being the intrinsic impedance of the free space. At 
resonance, the maximum stored electrical energy is the 
same as the maximum stored magnetic energy, and 
when the stored electrical energy is maximum, the 
stored magnetic energy is zero. Hence the stored energy 
in a section of an exponential line acting as a resonator 


can be obtained from its stored electrical energy alone. 
Stored electrical energy 


. lL pb(z) par 
<f if f | Inky. |°dzpdpdo 
2 0 a(z) 0 


l b(z) gaz sin? (BI — Bz 
Mex [ f 8 Be pdpdz 
0 a(z) ( ol 
p?{ In — 
a(z) 


M= 2) Bi. 


ll 


E,(p, 2) 


and 


E, = — 2iBetBle(a#!2) (22) 


U 


(23) 


where 


From the assumed variation of the characteristic imped- 
ance in an exponential line, 


Zo(z) = Zo(0)e%* 
(24) 
Substituting the results of (24) in (23), one obtains 


30erM? ( sin =) 


Z (0) 26 


U= (25) 


1 
RM? 3600 | Zo(0)d 1 ( 
bn? = Z?(0) 60 2a a? 
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For both series and parallel resonance, sin 26/=0, 


F 30er M?] 
Z(0) 
where € is the dielectric constant of the medium inside 


the transformer. Fig. 4 shows a comparison of the stored 
energy in uniform and exponential type resonators. 


U exponential res. / UV iniforms res. 
° SN 


x/B —> 


Fig. 4—Relative variation of stored energy in an 
exponential resonator. 


The wall losses consist of the losses resulting from the 
tangential component of the magnetic field at the lateral 
surface and at the end plates of the resonator. These can 
be determined from the surface integral of the tangen- 
tial H, over the entire surface. 

For parallel resonance, analogous to the \/4 type 
resonance in a uniform line resonator, the wall loss W 
becomes 


(26) 


x 
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where &, is the real part of the surface impedance of the 
metal forming the resonator. 
If Zo(0) be chosen 60 ohms and 


a—0 


then the Q of the resonator becomes 
2402? d oN 

Mine 
a b 


Rk, 
It is interesting to note that, when 


a—0 


the exponential-line resonator becomes a uniform-line 
resonator and the Q of such resonator® is 


2401? oN He 
(s t ae ) ; 
a b 


R, 
The characteristic impedance of the lines is assumed to 
be the same in both cases. 

The Q can be determined alternatively from the input 
reflection coefficient already derived. But the evaluation 
of QO from the field integrals reveals the characteristics of 
the resonator from the energy consideration. 

Fig. 5 shows a comparison of Q for uniform-line and 
exponential-line resonators for different values of a/B, 
when Z (0) =30 ohms. 


CONCLUSION 


An attempt has been made in this paper to describe 
a method by which the reflection coefficient of an ex- 
ponential-line transformer can be determined readily, 
from the Smith Chart, particularly when a frequency- 
sensitive load or a transmission line whose input im- 
pedance changes with frequency is terminated at the 


6 Ibid., p. 280. 
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Q 
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Fig. 5—Q gainin exponential resonator. Cavity length = — 


load end. Analyses of this nature may be helpful in ob- 
taining the optimum design parameters for the trans- 
former for any specific load. Also discussed is the possi- 
bility of using the exponential-line resonator to indicate 
how, for some range of Zo(0), the Q of an exponential 
resonator can be increased greatly in excess of what 
would be expected in a uniform-line resonator for the 
same type of resonance. Similar analyses can be made 
for other nonuniform line resonators.’ 

7 R.N. Ghose, “Synthesis of Nonuniform Line,” thesis submitted 


in partial fulfillment of requirements for degree of electrical engineer, 
Univ. of Ill.; 1956. 


Correction 


Tore N. Anderson, author of the paper “Rectangular 
and Ridge Waveguide,” which appeared on pages 201- 
209 of the October, 1956 issue of these TRANSACTIONS, 
regrets the omission of the following reference. The. 
illustration in Fig. 8 and the general equation for deflec- 
tion of pressurized waveguide were obtained from James 
L. Briggs and Joseph B. Brauer, Technical Note 
RADC-TN-54-10, p. 3; August, 1954. 
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Correspondence 


The Advantages of Expressing 
Standing-Wave Ratio in Decibels* 


It has been the practice at Wheeler Lab- 
oratories for the past ten years to express 
standing-wave ratio in decibels as opposed to 
a numerical ratio. So far as we are aware, 
this practice first evolved in the Bell Tele- 
phone Laboratories. The relationship is de- 
fined as follows: 


db SWR =20 logis VSWR. 


Expressing standing-wave ratio in db offers 
many advantages, as listed below. 

1) SWR may be read by simple sub- 
traction on a meter having a db scale. This 
avoids the necessity for precisely setting the 
level to unity on the ratio scale for the maxi- 
mum-voltage location on the standing wave. 

2) There is no need for designating SWR 


* Received by the PGMTT, March 14, 1957. 
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as a ratio of voltage (VSWR) or current or 
power, since decibels are so defined as to 
give the same number in any case. Similarly, 
expressing standing-wave ratio in db avoids 
the confusion between VSWR greater than 
unity, as is the custom in the U.S.A. and 
VSWR less than unity, as is the custom in 
Great Britain. 

3) When expressed in db, SWR goes to 
to zero for zero reflection. This feature is an 
aid in plotting and in thinking; it is espe- 
cially valuable in expressing tolerances. 

4) When expressed in db, SWR is pro- 
portional to magnitude of reflection coef- 
ficient within 4 per cent out to 6 db SWR 
(VSWR of 2:1). By comparison, a plot of 
VSWR vs reflection coefficient shows a de- 
parture from proportionality that is ten 
times greater at the same value. The closer 
proportionality of the db scale results in 
easier interpretation of data. 

5) The resultant of small reactive reflec- 


Sr 


July 


tions in a lossless line adding in the worst 
phase may be determined by merely adding 
individual values of SWR expressed in db. 
For example, if a line contains ten bumps, 
each having a reflection of 0.5 db SWR, the 
maximum possible reflection is 5.0 db SWR. 

6) Several computations involving re- 
flection coefficients are simplified when SWR 
is expressed in db. One such computation is 
the determination of minimum insertion loss 
by the loss-circle method. 

These advantages, borne out in our own 
experience developing microwave com- 
ponents and antennas, lead us to recom- 
mend a wider adoption of the practice of ex- 
pressing standing-wave ratio in db. 

David DETTINGER 
Wheeler Labs., Inc. 
Great Neck, N. Y. 


H. A. Wheeler and D. Dettinger, “Measuring the 
Efficiency of a Superheterodyne Converter by the In- 
put Impedance Circle Diagram,” Wheeler Monograph 
No. 9; 1949, 
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